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^ ! Abstract 



A new exact renormalization group equation for the effective average ac- 
tion of Euclidean quantum gravity is constructed. It is formulated in terms 
of the component fields appearing in the transverse-traceless decomposi- 
tion of the metric. It facilitates both the construction of an appropriate 
5^ \ infrared cutoff and the projection of the renormalization group flow onto a 

large class of truncated parameter spaces. The Einstein-Hilbert truncation 
is investigated in detail and the fixed point structure of the resulting flow 
is analyzed. Both a Gaussian and a non-Gaussian fixed point are found. If 
the non-Gaussian fixed point is present in the exact theory, quantum Ein- 
stein gravity is likely to be renormalizable at the nonperturbative level. In 
order to assess the reliability of the truncation a comprehensive analysis of 
the scheme dependence of universal quantities is performed. We find strong 
evidence supporting the hypothesis that 4-dimensional Einstein gravity is 
asymptotically safe, i.e. nonperturbatively renormalizable. The renormal- 
ization group improvement of the graviton propagator suggests a kind of 
dimensional reduction from 4 to 2 dimensions when spacetime is probed at 
sub-Planckian length scales. 



I. INTRODUCTION 



During the past decade, exact renormalization group equations in particular in the 
context of the effective average action , have become a powerful tool for the investigation 
of nonperturbative phenomena in both quantum field theory and in statistical physics. 
Those renormalization group (RG), or flow equations may be regarded the counterpart 
for the continuum of Wilson's renormalization group of iterated Kadanoff block spin 
transformations which had been formulated for discrete spin systems originally 0. In 
both cases the central idea is to "integrate out" all fluctuations with momenta larger 
than some cutoff k, and to take account of them by means of a modified dynamics for 
the remaining fluctuation modes with momenta smaller than k. This modified dynamics 
is governed by a scale dependent effective Hamiltonian or effective action, r^, whose 
fc-dependence is described by a functional differential equation, the exact RG equation. 

In quantum field theory this general strategy can be applied to both "effective" and 
"fundamental" theories. By definition, an effective theory is valid only if all relevant 
momenta of the process under consideration are close to some specific scale k which char- 
acterizes the theory. If is the action of an effective theory at scale k we can compute 
cross sections for the scattering of particles with momenta (or relevant momentum trans- 
fers) of the order of k, with all quantum effects included, by simply evaluating the tree 
diagrams of r^. Exact RG equations can be used in order to evolve to a smaller scale 
k' < k hj further "coarse graining" . 

Flow equations may also be used for a complete quantization of fundamental theories. 
If the latter has the classical action 5* one imposes the initial condition F^ = S* at the 
ultraviolet (UV) cutoff scale k, uses the exact RG equation to compute F^ for all k < k, and 
then sends k ^ and k ^ oo. Loosely speaking, the defining property of a fundamental 
theory is that the "continuum limit" k oo actually exists after the "renormalization" 
- in the traditional sense of the word - of finitely many parameters in the action; only a 
finite number of generalized couplings in Fq is undetermined and has to be taken from 
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the experiment. This is the case for perturbatively renormahzable theories [Q, but there 
are also perturbatively nonrenormalizable theories which admit a limit k ^ oo. The 
"continuum" limit of those nonperturbatively renormahzable theories is taken at a non- 
Gaussian fixed point of the RG fiow. It replaces the Gaussian fixed point which, at least 
implicitly, underlies the construction of perturbatively renormahzable theories Thus 
knowing its fixed point structure is crucial if one wants to assess whether a given model 
qualifies ClS db fundamental theory. 

In this paper we shall use a formulation where is the "effective average action" 0]. 
It is a coarse grained free energy functional which is constructed in close analogy with 
the standard effective action F to which it reduces in the limit of a vanishing infrared 
(IR) cutoff, k —>■ 0. The Euclidean functional integral for the generating functional W is 
modified by adding an IR cutoff term AkS to the classical action. It supplies a momentum 
dependent (mass)^-term TZkip^) for a mode of the quantum field with momentum p. 
The cutoff function TZkip^) vanishes for p"^ ^ fc^; hence the high-momentum modes get 
integrated out in the usual way. For p"^ <^ k"^ it behaves as TZkip"^) oc k"^ so that the small- 
momentum modes get suppressed in the path integral by a mass term oc A;^ . The scale 
dependent action F^ is closely related to the Legendre transform of the modified generating 
functional Wk- When regarded as a function of k, Tk runs along a RG trajectory in the 
space of all actions which starts at F-^ = 5* and ends at Fq = F. In the simplest case, the 
exact RG equation which describes this trajectory has the following structure: 



Here F^ denotes the infinite dimensional matrix of second functional derivatives of F^ 
with respect to all dynamical fields. 

This construction is fairly straightforward for matter fields, the inclusion of gauge 
fields introduces additional complications though. Using background gauge techniques. 



for gravity. Leaving the Faddeev-Popov ghosts aside, the effective average action for 




(1.1) 



a solution to this problem was given in refs. for Yang-Mills theory and in refs. 
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gravity, Tk[g^^,g^^], depends not only on the "ordinary" dynamical metric g^i, but also 
on a background metric g^i,. The conventional effective action r[g] is obtained as the 
k ^ limit of the functional Tk[g] = Tk[g,g = g] with the two metrics identified ^J^. 
The motivation for this construction is that in this manner Tk[g] becomes invariant under 
general coordinate transformations. 

Nonperturbative solutions to the above RG equation which do not require a small 
expansion parameter can be obtained by the method of "truncations". This means that 
one projects the RG flow from the infinite dimensional space of all actions onto some 
finite dimensional subspace which is particularly relevant for the problem at hand. In this 
manner the functional RG equation becomes an ordinary differential equation for a finite 
set of generalized couplings which serve as coordinates on this subspace. In ref. the RG 
flow of quantum General Relativity was projected on the 2-dimensional subspace spanned 
by the invariants / df^x^ and / d'^Xy/gR. This so-called Einstein-Hilbert truncation 
amounts to considering only functionals of the form 

^k[9,9] = (levrGfe)"^ J d'^x y/g [-R{g) + 2Afc| + classical gauge fixing. (1.2) 

Here Gk and are the running Newton constant and cosmological constant, respectively. 
More general and, therefore, more precise truncations would include higher powers of the 
curvature tensor as well as nonlocal terms [l^ which are not present classically. 

Quantum gravity is certainly a particularly interesting topic where exact RG equa- 
tions can lead to important new insights. As quantized Einstein gravity is perturbatively 
nonrenormalizable a natural option is to consider it an effective field theory [|lT|. Already 
within this setting quantum effects can be studied in a consistent and predictive way. In 
fact, in refs. ||r^ , p!3[] the running couplings Gk and obtained in were used to inves- 
tigate how quantum gravity effects modify the structure of black holes, and in [|l^ the 
implications for the cosmology of the Planck era in the very early Universe were studied. 
Along a different line of research it has been proposed |T5| that there are strong quantum 
gravitational effects also in the later stages of the cosmological evolution which even might 
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drive the cosmological constant to zero dynamically; the effective average action would 
be an ideal tool for exploring such infrared effects. 

An even more intriguing possibility is that, despite its perturbative nonrenormaliz- 
ability, quantized gravity exists nonperturbatively as a fundamental theory. It would then 
be mathematically consistent down to arbitrarily small length scales. A proposal along 
these lines is Weinberg's "asymptotic safety" scenario ITH]. It assumes that there exists a 
non-Gaussian RG fixed point at which the k oo limit can be taken, i.e. that the theory 
is "nonperturbatively renormalizable" in Wilson's sense. Asymptotic safety requires that 
the non-Gaussian fixed point is UV attractive (i.e. attractive for k oo) for finitely many 
parameters in the action, i.e. that its UV critical hypersurface is finite dimensional. This 
means that the RG trajectories along which the theory can flow as we send the cutoff k to 
infinity are labeled by only finitely many parameters. Therefore the theory is as predictive 
as any conventionally renormalizable theory; it is not plagued by the notorious increase 
of free parameters which is typical of effective theories. The set of generalized couplings 
for which the non-Gaussian fixed point is UV attractive should include the dimensionless 
Newton constant, g, and cosmological constant, A. 

Using the e-expansion, Weinberg showed already long ago that gravity in 2 + e dimen- 
sions (0 < £ ^ 1) is indeed asymptotically safe |T^. Further progress in this direction, in 
particular for d = 4, was hampered by the lack of an efficient calculational scheme which 
could be used to search for nonperturbative fixed points. 

As a solution to this problem which does not rely on the ^-expansion we propose 
to use the effective average action in order to find nontrivial fixed points {X*,g*,-'') 
of the gravitational RG flow. (The dots stand for the infinitely many other couplings 
which parametrize a generic action functional.) Using this approach, the case d = 2 + e 
was reanalyzed in a more general setting and, more importantly, it was shown that the 
Einstein-Hilbert truncation predicts the existence of a non-Gaussian fixed point {X*,g*) 
also in dimensions d > 2, in particular for d = 4 P JTB|JT7[] . 
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The crucial question which arises is whether this result is an artifact of the truncation 
used, or if it correctly reflects a property of the full theory. It is clear that in order to 
answer this question one would like to include further invariants into the truncation and 
to check whether the predictions stabilize. 

From the technical point of view such calculations are extremely complicated so that 
in the present paper we shall use a different method in order to get a first idea about 
the reliability of the nontrivial fixed point. We are going to analyze to what extent its 
location in the X-g plane and its attractivity properties (critical exponents) are scheme 
dependent. Here "scheme dependence" refers to the dependence on the cutoff operator 
TZk used in the derivation of the RG equation. 

First of all, TZk is a matrix in the space of irreducible component fields (see below) 
which is not uniquely determined by the general principles. Hence we can vary it to 
some extent. In fact, in the present paper we shall introduce a new cutoff ("cutoff of 
type B") whose matrix structure is different from the original one of ref. [^] ("cutoff of 
type A"). Either of these cutoffs is proportional to a "shape function" R^^^p^ /k'^) which 
describes the "thinning out" of degrees of freedom as we pass the threshold = k"^. Also 
this function can be varied in order to assess the scheme dependence of the fixed point 
properties. 

While in general only the critical exponents but not the location {X^,,g^, ■ ■ ■) of the 
fixed point are expected to be universal, i.e. scheme independent |T^, we shall argue that 
the product (7* A* is an observable quantity as well. For observables the 7?.fc-dependence is 
a pure truncation artifact; in an exact treatment all T^fc-dependencies cancel. The status 
of the Einstein-Hilbert truncation would be rather questionable if the fixed point was 
present for some cutoffs but absent for others. Instead, we find that it is actually there 
for all admissible cutoffs, and moreover that the observable g^,X^, is scheme independent 
with a quite unexpected precision. 

Our results strongly support the conjecture that the non-Gaussian fixed point is 
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present in the exact theory and is not a truncation artifact. Also another prerequisite of 
asymptotic safety turns out to be satisfied: we find that, for any cutoff, the fixed point is 
UV attractive in both directions of the X-g plane. 



Ultimately one would like to use more general truncations than ( |1.2|) in order to 
study the RG flow in a larger subspace. Typically this requires computations whose 

n 

algebraic complexity is quite formidable. Assume we make an ansatz Tj. = J2 

1=1 

containing n diffeomorphism invariant functionals Ii[g^u,gfiu]- In order to project the RG 
flow on the n-dimensional space with coordinates gj we must insert the ansatz into the 
RHS of the flow equation (|1 . 1|) . At this point the nontrivial problem, both conceptually 
and computationally, is to expand the trace with respect to a complete set of actions, 
{lj[gfiu,gfiu]}, in such a way that the Jj's retained in the ansatz are a subset of the j/s. 
The coefficients of the remaining j/s, those not present in F^, are set to zero by the 
truncation. In practice the projection on the {Jj}-subspace is done by inserting a set 
of metrics g^u^g^lu on both sides of eq. (|L1|) which give a nonzero value only to specific 
linear combinations of the Jj's. Provided one manages to compute the functional trace 
for sufficiently many g^^, ^^^.-pairs one can then deduce the ordinary differential equations 
for the generalized couplings giik). 

For the Einstein-Hilbert truncation this procedure is fairly simple since (ignoring the 
running of the gauge parameter) it is sufficient to insert for g^^, = g^y the metric of a 
family of spheres S'^ parametrized by their radius r. Their maximal symmetry facilitates 
the calculations considerably. With r kept as a free parameter, these metrics are general 
enough to disentangle / d'^x^/g oc r'^ and / d'^x^/gR oc r'^^^. But already when we in- 
clude invariants with four derivatives of the metric this method fails: the spheres cannot 
distinguish / d'^x^/gR'^ oc r'^~^ from / d'^x^/gR^^R^'^ oc r*^"^, for instance. 

These remarks hint at (at least) two major problems which one faces in general- 
izations of the exact RG approach to gravity, (i) The momentum dependent "mass" 
term Ai:S[hf^y,g^y] depends quadratically on the metric fluctuation h^^,, but also, via 
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IZk = 7^fc[(?] on the background g^^. In general it is a quite nontrivial task to construct 
a cutoff operator TZk[g] which has the desired properties mentioned above for a class of 
background metrics g^^ general enough for the projection on the truncation subspace. (ii) 
Assume we found an appropriate TZklg]- Then there arises the computational problem of 
evaluating the trace on the RHS of the RG equation for various fi'^zy's and ^^^'s. These 
metrics do not coincide when we allow for an evolution of the gauge fixing sector. Even if 
we ignore this complication, the Hessian under the trace is an extremely complicated 
nonminimal covariant matrix differential operator constructed from the curvature tensor 
and covariant derivatives D^. A priori, even for maximally symmetric backgrounds, not 
all derivatives are contracted to form powers of the covariant Laplacian = D^D^, 

(2) 

and is not diagonal in the space of fields with a definite helicity therefore. Hence 
standard heat kernel techniques or perhaps information about the spectrum of are of 
no help at this point. 

In this paper we outline a general strategy for tackling these problems. It is based 
upon York's "TT-decomposition" PD[ which is available on (almost) every spacetime 
manifold needed for our projection method. The idea is to decompose the fluctuation h^i, 
into a transverse, traceless tensor hj^^, a longitudinal-transverse tensor (parametrized by 
a transverse vector C,'^), a longitudinal-longitudinal tensor (parametrized by a scalar a), 
and a trace part (parametrized by another scalar (p). In the basis of the component fields 
{hj^j^, ^'^,a,(f)}, all -D^'s appear in powers of the Laplacian only, at least for the class of 
maximally symmetric backgrounds. The important point is that this decomposition can 

(2) 

be used in order to simplify the structure of Tl on essentially all backgrounds, not just 
on spheres. (On S"^ the TT-decomposition boils down to the familiar decomposition of 
h^i, with respect to pieces which are irreducible under the isometry group SO{d + 1). In 
some of the work following the original paper |Q this decomposition on S"^ had been used 
already [p!9|,pT|-p3|.) Compared to S'^ a certain complication arises, however, because the 
TT-decomposition is nonorthogonal in general. 
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The TT-decomposition also helps in solving the first problem, the construction of TZk, 
because A^S has a much simpler structure when expressed in terms of the component 
fields rather than the original h^^. 

This paper is organized as follows. 

In the first part (sections ^ and [TTID we describe the construction of a new RG equation 
where the component fields {/i^j,, a, 0} are used from the outset. Along the way we 
discuss the problems related to a proper identification of A^S. This part of the paper is 
meant to supply a set of tools which will become indispensable in future investigations 



when one includes further invariants into the truncation (i?^-terms for instance), if 
one adds matter fields, or if one allows for a running gauge fixing. 

As a first application, we revisit the Einstein-Hilbert truncation in the second part 
of the paper (sections and 0). We introduce a new cutoff TZk which is natural in 
the TT-language, and we use an arbitrary gauge parameter. This allows for a nontrivial 
comparison of the resulting RG equations and their fixed point properties with those of 
1^ whose cutoff operator TZk has a rather different structure. We find both a Gaussian 
and a non-Gaussian fixed point in the {g, A)-system and we perform a detailed analysis of 
their properties, in particular of their scheme dependence. The chances for realizing the 
asymptotic safety scenario in 4 dimensions will be discussed in detail. 

In section we investigate the implications of the non-Gaussian fixed point for the 
effective graviton propagator at large momenta. A kind of dimensional reduction from 4 
to 2 dimensions takes place in the vicinity of this fixed point. The asymptotic form of the 
propagator suggests that when 4-dimensional spacetime is probed by a very high-energetic 
graviton it appears to be effectively 2-dimensional. 

Various technical results, needed in the present paper but presented also with an eye 
towards future applications p4|, are relegated to a set of appendices. 

At this point the reader who is mostly interested in the results rather than their 
derivation can proceed directly to section |V|. 
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II. THE EXACT EVOLUTION EQUATION 



A. Gauge fixing 

Following IQ we define a scale dependent modification of the Euclidean functional 
integral for the generating functional by using the background gauge fixing technique 
QJ^. For this purpose we decompose the integration variable in the functional integral 
over all metrics, j^u, into a fixed background metric g^u and a fluctuation field h^j,^, 

= (jf^uix) + h^^{x) . (2.1) 

Then we replace the integration over 7^,^ by an integration over h^j^y. With the Faddeev- 
Popov ghosts and C*^ the generating functional may be written as 

Zfc [sources] = j Vh^^VC^VC^ 

y^ey.Y> -S[g + h]- Sgi[h]g]- SgY\h,C,C]g]- ^k^^ . (2.2) 

The first term in the exponential, 5* [7] = S[g + h], is the classical action which, for 
the moment, is assumed to be positive definite. It is invariant under arbitrary general 
coordinate transformations. S'gf denotes the gauge fixing term 

SAh-rg] = ^ j d'^x^r" F,[g,h]F,[g,h] . (2.3) 

It corresponds to the gauge condition F^[g, h] = 0. Linear gauge conditions, 

F,[g,h] = V2Kj^f[g]h^p, (2.4) 

are particularly convenient. In the present paper we use the harmonic gaugeQfor which 

^fm = Kr-'D, - \r'D, . (2.5) 



*For tfie flow equation in the conformal gauge (2D Liouville quantum gravity) see refs. [p^p 
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Here denotes the covariant derivative constructed from the background metric g^^, 
while we shall write for the covariant derivative involving the quantum metric 7^,^. In 
eq. (|2.4| ) we introduced the constant 



K = (327rG)-2 (2.6) 
where G denotes the bare Newton constant. The Faddeev-Popov operator associated with 



the gauge fixing (|2.4|) with ( |2.5D takes the form 

-M [7, = r'r'Dx {ip,D, + i.,D,) - r^r'DxiauD, . (2.7) 

It enters the functional integral (|2.2|) via the ghost action 



5gh[/i, C, C- -g] = -V2j d^x ^C,M[g + h, gf , C . (2.8) 

Furthermore, AkS and ^source are the cutoff and the source action, respectively. AkS 
provides an appropriate infrared cutoff for the integration variables and S'source introduces 
sources for the fields h^^, and C^. Their explicit structure will be discussed later on. 

B. Decomposition of the quantum fields 

For the calculations in the following sections it turns out to be convenient to decompose 
the gravitational field h^y according to (see e.g. [ ^ ) 

V = + ^M^^ + Dul^, + D^D^d - ^9f.uD^$ + ^g^y(j) . (2.9) 

To obtain this "TT-decomposition" one starts by splitting off the trace part /i^^ = g^u'P/d 
from h^j_y. It involves a scalar field (j). The remaining symmetric traceless tensor may be 
decomposed further into a transverse component K^^, and a longitudinal component h^^. 
Introducing a transverse vector field and another scalar a, the longitudinal tensor 
can be expressed by /i^^ = /iJJ + h^J^ with h^l = D^l^ + D^i^ and h^J^ = D^D^a - 
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g^yD'^a/d thereby ending up with eq. ( p.9|) . Thus the components of h^^ introduced by 
this transverse-traceless (TT-) decomposition obey the relations 

^^^/ij, = 0, D%, = 0, D^i, = 0, <P=-g,yh^''. (2.10) 

This decomposition is vahd for complete, closed Riemannian (i-spaces (i.e. compact 
Riemannian manifolds without boundary). As argued in PU|, its domain of validity 
can be extended to open, asymptotically flat ci-spaces, certain assumptions concerning 
the asymptotic behavior of the fields being made. From now on we assume that the 
gravitational background belongs to one of these classes of spaces. 

Obviously h^^ receives no contribution from those and a-modes which satisfy the 
Killing equation 



D^t + D,l, = Q (2.11) 



and the scalar equation 



D^D,a--^g^,D^a = {] , (2.12) 

respectively. Therefore such modes, referred to as unphysical and a-modes, have to 
be excluded from the functional integral. Considering the conformal Killing equation 

D^C, + DX^. - \g^.u DxC^ = (2.13) 

we recognize that the unphysical a-modes correspond to constants or are related via 
= D^d to proper conformal Killing vectors (PCKV's), i.e. solutions of eq. (|2.13|) 
which are not at the same time ordinary Killing vectors (KV's), |P7|| .p| 



tAs a consequence of the linearity of the conformal Killing equation we may add any Killing 
vector (which is always transversal) to its solutions and obtain another solution. For definite- 
ness we therefore define the PCKV's to be purely longitudinal. Then there is a one-to-one 
correspondence between the PCKV's and the nonconstant solutions of ( |2.12D . 
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By virtue of the decomposition ( p.9|) the inner product on the space of symmetric 
tensor fields may be decomposed according tof] 



cr 



(2) 



(2) 



:2.i4) 



From eq. (|2.14|) we see that, for a general background metric, only h^^, /i^^ and h^l 
form an orthogonal set, whereas and /i^^ are not orthogonal in general. This 
nonorthogonality manifests itself in the appearance of terms where the components 
and a mix. But at least for Einstein spaces, where R^i, = Cg^^p with C a constant, we 



find D^R^''' 



and therefore ( h , h 



,LT uLL\ 



AC j d'^x^a Df"^^ = 0. Thus 



{hj^i,, , h^^ , h^l} represents an orthogonal set of field components in this case. 

In order to determine the Jacobian Ji which appears in the functional integral (|2.2|) 
after performing the transformation of integration variables h^y — > {^^i/; ^' 
proceed as follows. We consider a Gaussian integral over /i^j, and reexpress it in terms of 
the component fields |27|: 



Vhny exp 



Here 



-2 (g^"'D^ + R^'^) -2R^'^Dx 

2DxR^'' ^{D'^f + DxR^PDp 



d'^x^ 




(2.15) 



(2.16) 



•^A remark concerning our notation: If not indicated otherwise each covariant derivative acts 
on everything that stands on the right of it. 
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is a Hermitian matrix differential operator. Since all functional integrals appearing in eq. 
( p.l5| ) are Gaussian they are easily evaluated. This leads to the Jacobian 

(2.17) 



= iVy'detfi^. 0) M^^'''^ ■ 

Here represents an infinite constant which may be absorbed into the normalization of 
the measure 'Dh^^. 

The notation adopted in eq. ( |2.17| ) has to be interpreted as follows. A prime at 
the determinant or the trace of an operator A indicates that all unphysical and a- 
eigenmodes of A, characterized by eqs. ( p.ll| ) and ( p.l2| ), are to be excluded from the 
calculation. A subscript at determinants or traces describes on which kind of field the 
operator A acts. We use the subscripts (0), (IT) and (2S'T^) for spin-0 fields a, transverse 
spin-1 fields and symmetric transverse traceless spin-2 fields hj^^^, repectively. The 
subscript (IT, 0) appearing in eq. ( ^.17] ) refers to a (c? + 1) x (d + l)-matrix differential 
operator whose first d columns act on traceless spin-1 fields C,'^ whereas the last column 
acts on spin-0 fields a. 

Likewise we decompose the ghost and the antighost into their orthogonal components 
according to 



:2.18) 



where Cj and C^^ are the transverse components of C„ and C^: D^Cj^ = 0, DnC^^ = 0. 



In order to compute the Jacobian J2 induced by the change of variables 
— > {C^^, fj} we write 



{Cj,r/} 



PC^PC. exp 



~ c,c 



= J2j VC^^" VCl Vf]Vfj exp - j (fx^ {CIC^^ + fi{-D^)ff] 

and perform the Grassmann functional integrals. The result is 

-1 



(0) 



(2.19) 



(2.20) 



In this case the constant f/-mode represents an unphysical mode which has to be excluded. 
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C. Momentum dependent redefinition of the component fields 



It will prove convenient to introduce new variables of integration, a, rj and rj, by 
means of the momentum dependent (nonlocal) redefinitions 



f] = ^-D'^f], 7] = ^-D^fj. (2.21) 

Here the operator Ric maps vectors onto vectors according to 

(mcvY = R'^'v^ . (2.22) 

Note that the transformations ( |2.21D are well defined and invertible since for any (physical) 
eigenmode the operators under the square roots of ( p.21|) have strictly positive eigenval- 



ues.0 This is due to the fact that these operators arise from the squares of /i^J and hj^^ 
and from (D^f/) (D^f^) by shifting all covariant derivatives to the right. Thus they cannot 
assume negative eigenvalues. For example, 



= ^ / d'^xVd^ ({D'y + -^D,R^^d\ a . (2.23) 



Furthermore, the spectra of the operators in ( p.21| ) do not even contain zeros, since the 
potential zero-modes coincide precisely with the aforementioned unphysical modes which 
have to be excluded. For instance, is a zero-mode of — Z)^ — Ric if and only if ^'^ is a 
Killing vector. 



^In order to make sure that the operators are indeed invertible we also assume that their 
eigenvalues do not have zero as an accumulation point. 
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Along the lines outlined in the previous subsection, we now determine the Jacobians 
for the transformation of integration variables. We obtain 



J. 



det'/ 1 



-D^ - Ric 



det;o) ( {D 



2\2 



a — 1 



= det/ 



(0) 



[2.24) 



for the transformations ^fj, ^ ^ and fj f], fj —>■ 7], respectively. (The integration 
measures have been chosen such that no additional infinite constants occur in the Jaco- 
bians.) The square brackets appearing in the subscript (1, [T]) at J3 indicate that the 
operator under consideration acts on spin-one fields which are transverse only for certain 
background metrics, because the property of transversality is not necessarily transmitted 
from C,'^ to However, at least for Einstein spaces is transverse as well. 

After carrying out this change of integration variables, Ji, J3 and J4 are the only 
Jacobians appearing in the generating functional since J2 and J5 cancel. 



D. The effective average action 

By adding an infrared (IR) cutoff A^S to the classical action under the path integral 
^ we obtain a scale-dependent generating functional Z^. The term A^S is chosen to 
depend on the fluctuation fields in such a way that their eigenmodes with respect to —D^ 
which correspond to large eigenvalues ^ k"^ are not influenced, whereas contributions 
from eigenmodes with small eigenvalues p"^ <^ k"^ are suppressed. In this sense Z^ de- 
scribes an effective theory at the scale k < k. For technical simplicity we implement the 
suppression of the low-momentum modes by momentum-dependent "mass" terms, i.e. by 
cutoffs which are quadratic in the fluctuation fields: 



A,S 



h, C, C;g]=^ {h, nr h) + {C, nf C) . (2.25) 
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Here the operators TZf^"^ and TZf^ are constructed from the covariant derivative with 
respect to the background metric, D^. Note that TZf'^'' and TZf^ must not depend on 
the quantum metric 7^,^ but only on the background metric g^u since otherwise the cutoff 
cannot be quadratic. In order to provide the desired behavior these operators must vanish 
for p^/k"^ — > 00 (in particular for k ^ 0) and must behave as TZk — ^ Zkk'^ for p'^/k'^ 0. 
(The meaning of the constant will be explained later.) As a consequence, all modes 
with <^ k'^ acquire a mass oc k. 

At this stage of the discussion it is not necessary to specify the explicit structure of 
the cutoff operators. We only mention the following point. According to appendix 
TZf'^'' and TZf" can be chosen such that, at the level of the component fields, 

A,s[h,C,C;g]=^ ^ (U,)^^^^ (2) + I E (V^i, (^^ ^^2) (2.26) 

with the index sets Ji = {/i^, ^, a, 0}, I2 = {C^ ,C,f],ri}. In contrast to generic cutoffs 
which are defined in terms of the component fields from the outset, the structure (|2.26|) 
allows us to return to the formulation in terms of the fundamental fields, eq. ( p. 251 ), in a 



straightforward way. (See appendix ^.) The set of operators ('^fc)f-^^2' ('^fc)v'iV'2 iii^^oduced 
by this realization of the cutoff may be fixed later on. Hermiticity demands that they 
satisfy (JZk)^^^^ = (J^k)^^^ and (7^fc)^^^^ = - (7^fc)J,^^2. Furthermore, (7^fc)^^^^ = if 
both ipi G {C^, T]} and 1^2 e {C^ , f]}, or if both ipi G {C^ , f]} and 1^2 e {C^ , f]}. 
A similar decomposition is applied to S'source- The source terms are defined as 

5'source[/^, C, C, J, K, K;g] = - (J, h) - (k, C) - {k, C) (2.27) 

with external sources J^^ , and for the fundamental fields C*^ and C^, re- 
spectively. Proceeding as described in appendix 0, an alternative form of S'source may be 
derived from eq. ( p. 271 ) where each component field is coupled to a certain component of 



the fundamental sources. Then, in terms of these "component sources , ^S'source takes the 
form 
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Ssouroeih, C, C, J, K,K;g]^-Y^ {J^, () - E 



Y: {K^,^) . (2.28) 



As a consequence, the functional 

Zfe [J, X, X; g] = Ji J3 J4 / Vhl, Vi, Va V(f> VC^^ PCj Vtj Vr] exp [ - g] 

-S^[h; g] - C, C; g] - A^S [h, C, C; g] - S^^^^h, C, C, J, K, K; g]^ (2.29) 

as well as the scale-dependent generating functional for the connected Green's functions, 

(2.30) 



wJj,K,K;g] =lnZ, 



J,K,K;9 



may be viewed as functional of either the fundamental or the component sources. Fur- 
thermore, we may derive /c-dependent classical fields for both fundamental and component 
fields in terms of functional derivatives of Wk- In either case the fc-dependent classical 
fields represent expectation values (q) of quantum fields q, in the sense that all degrees of 
freedom corresponding to momenta with > /c^ have been averaged out. The classical 
fundamental fields are given by 



1 SWk 



C\- 1 ^^'^ v^^iC^)- ^ 



(2.31) 



(2.32) 



and the classical component fields are obtained as 

" = ■ 

Here wc arc making use of the shorthand notation x = {h^ ^^^'^^'P^C'^ -.C^ ^V^V) the 
quantum component fields, J = (J/jT, J^, J^., J^, Kqt, Kqt, Kfj, K^^) for their sources and 
(/? = iji^ , ^, CT, 0, , v"^, g, g) for the classical component fields. We may reconstruct the 
classical fundamental fields from </? according to 



h 



-D^ - Ric 



-D^ - RicJ ' + -i9^i.v(t> 



-0 - \9,.& 



d 



g . 



(2.33) 
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Performing a Legendre-transformation on Wk with respect to J^,y, K'^ and leads 
to the following scale-dependent modification of the effective action: 

(2.34) 



Tk h,v,v;g = Uh) + (K,v) + {K,v) - Wk J,K,K;g 



Since eqs. (p^), {^^) imply 



{J,^) = Uh) + (K,v) + {K,v) 



(2.35) 



it is clear that eq. (|2.34| ) also would result from Legendre-transforming Wk with respect 
to the component sources. Denoting the corresponding Legendre-transform by r^°™'' we 
have T'f^'^^lh^ , C,, a, (p, v'^, v'^, g, g; g] = Tk[h, v, v; g] where the arguments of F™™'' and Ffc 
are related by (|2.33|) . 

The effective average action proper, Tk, is defined as the difference between f ^ and 
the cutoff action with the classical fields inserted P^J5|]: 



Tk[g, g, V, v] = - g, v, v; g] - AkS[g - g, v, v; g] . 



(2.36) 



Here we expressed h^i, in terms of the classical counterpart g^^ of the quantum metric 
1fj.u = g^lu + h^y which, by definition, is given by 



^2.37) 



The main advantage of the background gauge is that it makes Tk a gauge invariant 
functional of its agruments [Q. It is invariant under general coordinate transformations 
of the form 



rfc[$] = rfc[$ + , $ = {g^y, g^y, v", V, 



(2.38) 



where Cu is the Lie derivative with respect to the generating vector field u^{x). Since 
general coordinate invariance ensures that no symmetry violating terms occur in the 
course of the evolution of Tk the class of consistent truncations is restricted to those which 
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involve only invariant field combinations. This is important for practical applications of 
the evolution equation. 

We are mainly interested in the exclusively (y'^i.-dependent functional 



TfcM = ^k[9, 9,0,0] 



(2.39) 



In the limit — > it coincides with the conventional effective action r[(7^j,], the generator 

of the IPI graviton Green's functions M: r[g] = limTi^lg]. However, in order to derive 

fc— >o 

an exact evolution equation it is necessary to retain the dependence on the ghost fields 
and g^^. 

From the definition of the effective average action it follows that satisfies the integro- 
differential equation 

exp {-Tkig, g,v,v]} = J ^^^7^ ^^C, exp [ - S[h, C, C; g] 



5v,. 



with 



X exp |— AfcS* h — g + g,C — v,C — v; g | 



S[h, C, C- -g] ^ S[g + h] + S^,[h- g] + S^^[K C, C- g] . 



(2.40) 



(2.41) 



Eq. ( p.40| ) may be derived by inserting the definition of F^ into ( p.2|) and replacing the 
sources according to 

1 Sfk 1 6fk 1 5ffc 



[2A2) 



E. Derivation of the exact evolution equation 

The exact renormalization group equation describes the change of the action functional 
Ffc induced by a change in the scale k. It may be obtained as follows. Differentiating the 
functional integral (|2.29|) with respect to t = In A; leads to 
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d^Wk = -Tr' 
2 



Cl,C2 6-fl 



+ -Tr' 
2 



^2.43) 



Here we used eq. ( p.30| ) and adopted the matrix notation on the RHS of ( |2.43| ) which, in 
turn, can be expressed in terms of the Hessian 

1 



[il. 



^2.44) 



/g{x)g{y) ^^i{x)5^j{y) 

with [j] = for commuting fields ipj and [j] = 1 for Grassmann fields ipj. Since the 
connected two point function 

1 



(G'fc)ij(x, y) = (xiix) Xj{y)) - <^i(a;) '^j{y) 



and r^^^ are inverse matrices in the sense that 



(2.45) 



/ A ^) {GkUx, y) {ftYiy, A = S„ 



I S{ x-z) 



(2.46) 



we may replace the expectation values {Xii^)Xjiy)) appearing in eq. ( |2.43|) with 
(r^ )~- {x,y) + ipi{x)ipj{y). Then performing a Legendre-transformation according to 



eq. ( |2.34] ) and subtracting the cutoff action AkS[h,v,v; g] yields the desired exact renor- 
malization group equation: 



dtTk [g,g,v,v] = -Tr' 



1. 



.Cl,C2 6/l 



+ -Tr' 
2 



E {Tf\grg,v,v] + nky'd,ink)^,^^^ 



(2.47) 



Here we wrote (7^fc)cif2 = ('^fc)(Ci>(C2>> (^fc)v'iV2 = (^fc)(V'i>(V'2> and introduced the index 
sets 



h = {h ,^,cr,0}, l2 = {v ,v ,g,g} . 



(2.48) 



In a position space representation, the operators appearing on the RHS of the flow 
equation are given by matrix elements whose traces are evaluated according to 
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for instance. The notation adopted for the matrix elements is similar to eq. (|2.44|) ; for 
example, 

By virtue of the properties of TZk discussed above the traces appearing in the flow equation 
( p.4?| ) are perfectly convergent for all values of k < k. 

Provided we impose the correct initial condition at the UV scale k = k we can, in 
principle, determine the functional integral ( p.2|) by integrating the flow equation from 
k down to k and letting k ^ 0, k oo after appropriate renormalizations. The initial 
condition F-^ can be obtained from the integro-differential equation ( ^^.40 ). For sufficiently 



large values of k, the cutoff term in eq. ( |2.40D strongly suppresses fluctuations with 



{h, C, C) 7^ {h, V, v) so that the main contribution to the functional integral results from 
small fluctuations about {h,C,C) = {h,v,v). This field configuration corresponds to the 
global minimum of the total action in the exponential of eq. (|2.40| ). Performing a saddle 



point expansion of the functional integral about this minimum leads to 

Tk[g, -g, V, v] = S[g - g, v, v;g]-^ In Det' [S^'^ + n^) (2.51) 

where the second term contains one-loop effects. For k = k oo they amount to an often 
unimportant shift in the bare parameters of S which can be ignored usually. For finite 
k, additional contributions from the determinant occur which are suppressed by inverse 



powers of k [^. Therefore we obtain the initial value for /c — > oo 

^[9, 9, V, v] = S[g] + Sg{[g -g;g] + 5gh[^ - 9, v, v; g] . (2.52) 
At the level of the functional rfc[5f^i,] this initial condition boils down to 

h[9] = S[9] . (2.53) 
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So far we assumed the fundamental action to be positive definite. However, the 
Einstein-Hilbert action, for instance, does not have this property which is due to the 
appearance of a "wrong-sign" kinetic term associated with the conformal factor. In such 
cases it is nevertheless possible to formulate a well-defined evolution equation if the signs 
of the cutoff operators TZk are properly adjusted 0. We will return to this point in the 
next section. 



F. A special case: Einstein backgrounds 

Before continuing we summarize the simplifications that occur for Einstein back- 
grounds, for which Rf^i, = Cg^p with C a constant. In this case the decomposition 
( p.9|) of h^i, is completely orthogonal. In fact, thanks to the Einstein condition, the ^^-a 
mixing terms in the inner product ( |2.14] ) vanish so that K^^,, /i^J, h^^ and hj^l form an 
orthogonal set. 

Furthermore, for Einstein spaces the Jacobians appearing in the path integral ( |2.29|) 
cancel, at least up to an (infinite) constant which can be absorbed into the normalization 
of the integration measure. This can be seen as follows: 



Ji = iV^det(i2._o) M^^'"^ 

J Vl, Va exp [ - I rf^x { - 2^, {d' + C) C + a (^(D')' + C dA a^J^ 



iVi^/det'd^) i^-D^-C)Jdet[,^i-D 



)2 



^ det'(o) - = N, J,' (2.54) 



Here A^^i and N2 are unimportant constants so that Ji J3J4 is indeed field independent. 

Finally, for Einstein spaces the field redefinitions in the gravitational sector take the 
form 



r = \/-D' -Ci^, a = 4^N-D^-^^ . (2.55) 
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As a consequence, we find that D^^^^ = \J —D'^ — 2C D^^^^. Thus, transversahty of 
imphes that is transverse as well. 

III. TRUNCATIONS AND CUTOFFS 

A. A general class of truncations 

In practical applications of the exact evolution equation one encounters the problem 
of dealing with an infinite system of coupled differential equations since the evolution 
equation describes trajectories in an infinite dimensional space of action functionals. In 
general it is impossible to find an exact solution so that we are forced to rely on ap- 
proximations. A powerful nonperturbative approximation scheme is the truncation of the 
parameter space, i.e. only a finite number of couplings is considered. In this manner the 
renormalization group flow of Tk is projected onto a flnite-dimensional subspace of action 
functionals. In practice one makes an ansatz for that comprises only a few couplings 
and inserts it on both sides of eq. ( |2.47D , thereby obtaining a truncated evolution equa- 
tion. By projecting the RHS of this equation onto the space of operators appearing on 
the LHS one arrives at a set of coupled differential equations for the couplings taken into 
account. 

As discussed in refs. [^^, Ward identities provide an important tool for judging the 



admissability and quantitative reliability of a given truncation; approximate solutions of 
the flow equation are not necessarily consistent with the Ward identities, in contrast to 
the exact solution. Therefore, only those truncations which are indeed consistent with 
the Ward identities, at least up to a certain degree of accuracy, will yield reliable results. 
The Ward identities to be considered here are modifled by additional terms coming from 
the cutoff which are not present in the ordinary identities. Since AfcS* vanishes as /c — 
the ordinary Ward identities are recovered in this limit. 

In ref. H the modifled Ward identities were derived for the gravitational effective 
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average action Tk[g, g, v, v; (3, t] where f3 and r are auxiliary sources for the BRS variations 
of the graviton and the ghosts which are needed in order to formulate the Ward identities. 
Setting /5 = r = in the argument of this more general functional we obtain the action 
Tk[g,g,v,v] discussed in the present paper. (It would be straightforward to include the 
/3,r-sources also in the new formulation of the flow equation, but we will not need them 
in the following.) 

In the Ward identities were used to test the consistency of truncations of the form 
Tk[g,g,v,v] = Tk[g] + fk[g,g] + Sgi[g -g;g] + Sgh[g - g,v,v;g] (3.1) 



with Tk[g] defined as in eq. (|2.39| ). The term ric[g,g] encodes the quantum corrections 



to the gauge fixing term. This interpretation of rfc[(7,(7] is obvious because for g ^ g 
the purely gravitational part of eq. ( |3.1D implies Tk[g,g, 0,0] — Tk[g, g, 0,0] = Tk[g,g] + 
Sgi[9 ^ QiO]- By definition, Ti:[g,g] = 0. In the ansatz (|3.1|) the ghost dependence has 
been extracted in terms of the classical Sg^, thereby neglecting the evolution of the ghost 
action. This guarantees that the initial condition ( |2.52| ) is satisfied automatically in the 
ghost sector. In the gravitational sector it requires F-^ = S", F-^ = 0. For truncations of 
the type ( p.l| ) the Ward identities demand that Tk[g] is a gauge invariant functional of g^^ 
and they yield a constraint equation for rk[g,g]. To lowest order, this equation is solved 
by Ffc = < k. In the Einstein-Hilbert truncation we go beyond this approximation 
and set F^. oc Sgf with a constant of proportionality which vanishes aX k = k; it takes the 
running of the graviton's wave function normalization into account (see below). 

Inserting the ansatz (|3.1|) into the exact evolution equation ( p.47|) leads to a truncated 
renormalization group equation which describes the evolution of F^ in the subspace of 
action functional spanned by (|3.1D. The equation governing the evolution of the purely 
gravitational action 

^k[9,9] = Ffc[^,^,0,0] =fk[g] + Sgi[g-g]g]+fk[g,g] (3.2) 

reads 
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^ )'llj2lpl 



(3.3) 



Here r^.^-* and S^^ are the Hessians of Tk[g,g] and Sgii[h,v,v, g] with respect to the grav- 
itational and the ghost component fields, respectively. They are taken at fixed g^^. 



:(2) 



B. Specification of the cutoff 

In order to obtain a tractable evolution equation for a given truncation it is convenient 
to use a cutoff which is adapted to this truncation but still has the general suppression 
properties described in subsection PD| . It is desirable to start from a definition of A^^* that 



brings about the correct suppression of low-momentum modes for a class of truncations 
and of gravitational backgrounds which is as large as possible. 

A convenient, adapted cutoff can be found by the following rule [§,0. Given a 



truncation, we assume that for g = g the kinetic operators of all modes with a definite 
helicity are of the form {T^^^)ij = fij{—D^, k, . . .) where {fij} is a set of c-number functions 

(2) 

and the indices i,j refer to the different types of fields. (The difficulty of bringing 
to this form is one of the main reasons for using the TT-decomposition. At least for 
maximally symmetric spaces it allows us to eliminate all covariant derivatives which do 
not appear as a Laplacian = g^'^D^Du.) Then we choose the cutoff in such a way that 
the structure 

(rf + 7^,),, = [-D^ + eR^'\-D^/e), k,...) (3.4) 

is achieved. Here the function R^'^\y), y = —D^/k'^, describes the details of the mode sup- 
pression; it is required to satisfy the boundary conditions R^'^^{0) = 1 and lim R^'^^y) = 0, 



but is arbitrary otherwise. By virtue of eq. ( |3.4| ), the inverse propagator of a field mode 
with covariant momentum square = —D^ is given by -|- k"^ R^^^p"^ / k"^) which equals 
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for ^ k'^ and + k'^ for ^ A;^. This means that the small- J9^ modes, and only 
those, have acquired a mass oc k which leads to the desired suppression. 

In the next section we shall see in detail that for the truncations used in the present 
paper we can comply with the above rule by using the following cutoff operator: 



Zi^K^r^k^R^'^ (-D'/k 



d 



d-1 



Here Pu is defined as 



P^ = ^D^ + k^R^'\-Dye) . 



(3.5) 



(3.6) 



The remaining cutoff operators not listed in eq. ( |3.5| ) are set to zero. The Z^s are 
constants which, again by using ( p.4|) , will be fixed in terms of the generalized couplings 
appearing in the ansatz for F^. The cutoff ( p.5|) is inspired by the TZk used in |^ for S"*^. 

If eq. ( ^.41 ) allows us to choose Z^'^ > for all ( G , ^ , a , (f)} and Zf'^ = one 
obtains a positive definite AkS in the gravitational sector. In this case exp(— A^S*) is a 
damped exponential which indeed suppresses the contributions from the low-momentum 
modes. In the following sections we shall focus on the Einstein-Hilbert truncation for 
Tk which suffers from the conformal factor problem: its kinetic term for (f) is negative 
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definite. As a consequence, eq. ( ^.41 ) forces us to work witli a < 0. Hence, in tlie 
0-sector, AkS is negative definite and, at least at a naive level, exp(— A^S") seems to 
enhance rather than suppress the low-momentum modes. As we discussed in detail in ref. 
0] we nevertheless believe that the rule (|3.4| ), i.e. allowing for Z^'^ < 0, is correct also 
in this case. We emphasize that the RHS of the flow equation, contrary to the Euclidean 
path integral, is perfectly well-defined even if 5* and are not positive definite. 

At this point it should be mentioned that the situation with respect to the positivity 
of the action improves considerably by including higher- derivative terms in S and the 
truncated Tk since these actions are bounded below, provided we choose the correct sign 
in front of these higher-derivative terms. Furthermore, their quadratic forms are positive 
definite at least for sufficiently large momenta, and so is the cutoff. For a study of the 
evolution equation for i?^-gravity we refer to |2|]. 

As compared to the original paper |^, the cutoff ( |3.5| ) has a rather different structure 
which is due to the fact that it is formulated in terms of the component fields arising 
form the TT-decomposition. Contrary to the original one of ref. 0, the new cutoff (|3.5|) 
is defined for all values of a. This is one of the main advantages of the new approach. 

Note that in refs. [22,1£] where the TT-decomposition was used on S'^ the actual 
construction of the effective average action and its RG equation was omitted and has 
been replaced by an ad hoc modification of the standard one-loop determinants. No 
AfcS* has been specified at the component field level. Hence the scale dependent action 
constructed in this manner has no reason to respect the general properties of an effective 
average action 0. Despite the use of the component fields in [ p2| , p!^ their cutoff seems 
to be more similar to the original one in than to the new one of the present paper. In 
fact, it represents an a-dependent generalization of the cutoff in 0], in the sense that the 
latter is recovered from the one of refs. |2^,TD] by setting a = 1. 

From now on we will refer to the cutoff used in the original paper |^ and in 
as the cutoff of type A. However, one has to keep in mind that the existence of a 
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corresponding A^S is guaranteed only for a = 1, i.e. the case considered in |]6|Jlj 
Furthermore the cutoff ( p.5| ) of the present paper will be referred to as the cutoff type 
B; it is defined for all values of a. 

Each cutoff type contains the shape function R^^\ A particularly suitable choice is 
the exponential shape function 



R^'\y)=y[exp{y)-1]-' . 



(3.7) 



In order to check the scheme independence of universal quantities we employ a one- 
parameter "deformation" of ( p. 7] ), the class of exponential shape functions, 



s) = sy [exp{sy) - 1] \ 



with s parametrizing the profile of R^^^ |T9|. Another admissible choice is the following 
class of shape functions with compact support: 



R^'\y-b) 



exp 




{y- 1.5) ^exp (b-y) 



y<b 

b <y <1.5 
y>1.5 



(3.9) 



Here 6 G [0, 1.5) parametrizes the profile of R^^\ 

For our analysis of the flow equation in section we shall use both cutoff types with 
both classes of shape functions. 



IV. THE EINSTEIN-HILBERT TRUNCATION 



A. The ansatz 



In this section we use a simple truncation to derive the renormalization group flow 
of the Newton and the cosmological "constant" by means of the truncated flow equation 
( ^.31 ). In our example we assume that, at the UV scale k oo, gravity is described by 
the classical Einstein-Hilbert action in d dimensions. 
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hig] = S[9] = J d''x^{-R{g) + 2X} . (4.1) 

For the investigation of the evolution of rk[g,g] towards smaller scales k < k we consider 
a truncated action functional of the following form: 

rk[g,g] = 2K^ZNk J d''x^{-Rig)+2Xk} 

J d'x ^gr" {:Ffg^,) (^r^7p.) • (4-2) 

Eq. (^4.2|) is obtained from S + S^i by replacing 

G Gk = Z^l G , A Afc , a Z^la (4.3) 



so that its form agrees with that of the gravitational sector of the ansatz ( p.l|) with 

fk[g, g] = '^'^^^^ / d'x Vdr-" (-^r^p.) • (4.4) 

Generally speaking also the gauge fixing parameter a should be treated as a running 
quantity, a a^. Fortunately there is a simple shortcut which avoids an explicit com- 
putation of the corresponding /3-function. In fact, there are general arguments showing 
that a should have a (IR attractive) fixed point at a* = 0. This means that the initial 
condition = leads to = for all k < k. Thus, even using the truncation with a 
constant a, we can take the correct "flow" of the gauge fixing term into account simply 
by setting a = 0. 

In Yang-Mills theory the existence of the fixed point a^, = has been demonstrated 
for a truncation containing a covariant gauge fixing [2£], while for the axial gauge a 
nonperturbative proof is available [0 . The following general argument^ suggests that 
this fixed point should exist in any gauge theory, including gravity In the ordinary 



functional integral, the limit a — > corresponds to a sharp implementation of the gauge 



**We are grateful to J. M. Pawlowski for a discussion of this point. 
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fixing condition, i.e. exp(— S'gf) becomes proportional to S[F^]. The domain of the / Vh^^- 
integration consists of those /i^^'s which satisfy the gauge condition exactly, = 0. 
Adding the IR cutoff at k amounts to suppressing some of the /i^,y-modes while retaining 
the others. But since all of them satisfy = 0, it is clear that a variation of k cannot 
change the domain of the /i^i^-integration. The delta-functional S[Ffj] continues to be 
present for any value of k if it was there originally. Hence a vanishes for arbitrary k. 



B. Projecting the flow equation 



The fc-dependent couplings in eq. ( ^4.3| ) satisfy the initial conditions = A and 
Z^-j: = 1 which implies = G. Here the UV scale k is taken to be large but finite. The 
evolution of Zpfk and A^ towards smaller scales may now be determined as follows. As a 
first step the ansatz ( [4. 21 ) is inserted into both sides of the truncated flow equation (^3.31). 
Then we may set g^i, = g^^. As a consequence, the gauge fixing term drops out from the 
LHS which then reads 

dtTk[9, g] = 2k^ I d^x ^ [-R{g) dtZ^k + 2dt {z^^Xk)] . (4.5) 

Performing a derivative expansion on the RHS we may extract those contributions which 
are proportional to the operators spanning the LHS, i.e. J d'^x^ and J d'^x^R. Then, 
comparing the coefficients of these operators yields a system of coupled differential equa- 
tions for Zj\fic and A^. It describes the projection of the renormalization group flow onto 
the two-dimensional subspace of the space of all action functionals which is spanned by 
/ d^Xy/g and / d'^x^R. 

It is important to note that during this calculation we may insert any metric g^v 
that is general enough to allow for a unique identification of the operators / d'^x^/g and 
/ d'^x^R. In practice it proves particularly convenient to exploit this freedom by choosing 
the gravitional background to be maximally symmetric. Such spaces form a special class 
of Einstein spaces and are characterized by 
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R R 

with the curvature scalar R considered a constant number rather than a functional of the 
metric. In the following we restrict our considerations to maximally symmetric spaces 
with positive curvature scalar ^ > 0, i.e. (i-spheres S'^. For d fixed, S'^ is parametrized 
by the radius r of the spheres, which is related to the curvature scalar and the volume in 
the usual way, 

R=^^^-^, [d'xV^=^iAnrY^ . (4.7) 



r2 '7 T{d) 

Before continuing with the evaluation of the RHS of the flow equation we have to 
comment on the properties of fields defined on spherical backgrounds. According to 
appendix 0, we may expand the quantum and classical component fields in terms of 
spherical harmonics Tj^, T^™, T^™, which form complete sets of orthogonal eigenfunctions 
with respect to the corresponding covariant Laplacians. The expansions of hj^^, 0, C^, 



and their classical counterparts can be inferred directly from eq. (^^). The remaining 
component fields are expanded according to 

oo Dt(d,l) oo Dt(d,0) 

1=2 m=l 1=2 m=l 

oo Di{dfi) oo Diidfi) 

vix)=Y: E VimT'^'^ix), r^(x)=E E Virr.T'^'ix). (4.8) 

1=1 m=l 1=1 m=l 

Similar expansions hold for the associated classical fields. 

Note that in eq. ( ^.8| ) the summations do not start at / = 1 for vectors and at / = for 
scalars as in eq. ( p3D , but at / = 2 for and a, and at / = 1 for the scalar ghost fields. 
The modes omitted here are the KV's (T^^^'™'), the solutions of the scalar equation (|2.12|) 
which are in one-to-one correspondence with the PCKV's (T'^^'™), and the constants 
^rpi=o,m=iy As discussed in subsection [11 B| , the fundamental fields obtain no contribution 
from these modes. Therefore they have to be excluded from eq. 
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This exclusion is also of importance for the momentum dependent field redefinitions 
( p.21| ) because they would not be well-defined otherwise, as can be seen e.g. from 

oo Dt(d,0) 

a{x) = E E ; ^'"^(^) • (4.9) 

1=2 m=i yA^(d,0)(A,(rf,0)-^) 

Eq. (19) follows from inverting eq. ( 2.21 ) and then inserting eq. (|4.8|) . The eigenvalues 
corresponding to the modes excluded, i.e. Ao((i, 0) = and Ai((i, 0) = (see table 1 
in appendix 0), would lead to a vanishing denominator in eq. ( [4.9| ). Similar arguments 



hold for the other fields in ( |2.21 



We may now split the quantum field (f) into a part (pi spanned by the same set of 
eigenfunctions as a, and a part (pQ containing the contributions from the remaining modes: 

1 Diidfi) oo Di(d,0) 

0(x) = 0o(x) + 0i(x), 0o(a;) = E E 0/mT'™(x) , 0i(x) = E E </'«-T'™(x) (4.10) 

1=0 m=l 1=2 m=l 

The orthogonality of the spherical harmonics implies (01, 0o) = {<^,4>o) = so that 
{4>,(f)) = (00) 0o) + (01) 0i) and in particular (a, 0) = (cr, 0i). As a consequence, decom- 
posing according to eq. ( [4.10|) ensures that any nonzero term bilinear in the scalar fields 
is of such a form that the scalars involved are spanned by the same set of eigenfunctions. 
The same is true for the corresponding classical fields 0o and 0i. 



C. Evaluation of the functional trace 

Let us now return to the evaluation of the flow equation. On the RHS we need the 
operator F^^-* [g, g]. For our purposes it is sufficient to determine this operator at g^^i, = g^^. 
It may be derived by expanding Tk[g,g] according to 

F,[^ + h,g]= T,[-g,g] + OQi) + TT'^[h--g] + 0{¥) (4.11) 

and retaining only the part quadratic in li^y^ i.e. T^^\h]g]. For our truncation it takes 
the form 
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pquad 
k 



h;g= K^Z^k / d'^x ^ h 



1 10 



1 



1 — a r 



a 



(4.12) 



where g^^ is fixed but still arbitrary. In order to (partially) diagonalize this quadratic 
form we insert the family of spherical background metrics into eq. ( [4.12|) and decompose 
h^u according to eq. ( p.33| ). Then we use the classical analog of eq. ( |4.10| ) to decompose 
as well. This leads to 



+ -^^ -D'^ + Av{d,a)R-2a\k + Cs2{d,a){a -D^ + As2{d, a)R + Bs2id, a)\ 



a 



a 



+2Cs3id,a) 



R 



d-1 



-a 



4'&{4'o,<t>i} 



-D^ + Asi{d, a)R + Bsi{d, a)X^ 



(4.13) 



Here the A's, -B's, and C"s are functions of the dimensionality d and the gauge parameter 
a. The explicit expressions for these coefficients can be found in appendix 0. 

Note that this partial diagonalization simplifies further calculations considerably, and 
this is the main reason for using the decomposition ( ^.33p and specifying a concrete 
background. In contrast to the case a = 1 considered in a complete diagonalization 
cannot be achieved by merely splitting off the trace part from h^^, since eq. ( 4.12|) contains 
additional terms proportional to 1 — a which introduce mixings between the traceless part 
and (j). To be more precise, it is the term / d'^Xy/gh^i^[g'^'^ D^D^ — 6^3" Dp]hP'^ that gives 
rise to such cross terms. 

In terms of the component fields these cross terms boil down to a purely scalar (x-c^ 
mixing term that vanishes for the spherical harmonics 2"'=o.™=i and T'"^'™. Since these 
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modes contribute to (but not to a) we cannot directly invert the associated matrix 
differential operator ((^k^^^)ij^ . . _ _ _ . As a way out, we split according to eq. ( |4.10| ) 



into 00 and 0i. This has the effect that only mixings between the scalars a and 0i survive, 
which have the same set of eigenfunctions T'™" starting at / = 2. Hence the resulting matrix 
differential operator ({T^?)i^^ - - - is invertible, but since this additional split of 

affects the matrix structure of this operator it leads to a slightly modified flow equation. 



In fact, on the RHS of eq. (|3.3| ) the summation in the gravitational sector now runs over 
the set of fields {W ,^,(f)o,a,(f)i}, with {Jlkj^o^o = (^fc)</.i<^i = (T^fe)^,/.- 

In the context of the Einstein-Hilbert truncation it is only the a-dependence that in- 
troduces mixings of and the traceless part of h^u and therefore necessitates the decom- 
positions ( |2.9| ), ( |4.1(]|) . In general the inclusion of higher derivative terms like / d'^x y/gR^ 
and of matter fields leads to similar mixings. 

In order to determine the contributions from the ghosts appearing on the RHS of eq. 
( p.3|) we set g^i, = gfj,u in S'gh and assume that gfj,u corresponds to a spherical background. 
Then we decompose the ghost fields according to eq. ( |2.33| ) which leads to 



S'gh [0, v,v;g] = V2 J d'^x ^/g 



d. 



-D^ - 2- 

d. 



(4.14) 



From now on the bars are omitted form the metric, the curvature and the operators 
and Pfc. Note that the decomposition of the ghosts is not really necessary, but it allows 
for a comparison with the results obtained in 

(2) 

At this point we can continue with the adaption of the cutoff to the operators T), and 



of eqs. ( [4.131 ), (|4.14| ). According to the rule (|3.4| ) the Z^s have to be chosen as 



Zp" = Z^fc, zi^=-ZNk, Zt"" = Cs2id,a)Cs3{d,a)ZNk, Z^ = Cs2{d,a)Z^ 



a 



Nk 



Zt^' = Zt^^ = Cs2id, a)CsM o^)ZNk , Zf^^ = Zf = V2 . (4.15) 

— (2) 

Thus, for g = g the nonvanishing entries of the matrix differential operators + TZ^ 



and Si^^ + TZk take the form 
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[Tt\g.g]+n, 
Vf[g,g]+n 



hTJiT 



55 



ZNkfi'^ — Pk + Av(d,a)R- 2aXk 
a L 



ZNkti^ Cs2{d, a) Pk + As2{d, a)R + Bs2{d, a)Xi, 



Tf[g,g] + n, 



o"<pi 



ZNki^^ Cs2{d, a)Cs3{d, a) J Pk\ Pk 



R 



d-l 



{T^^\g,g] + nu 



</'0<^0 



(rf|j,j| + R.)._._ 

^ATfc/t^ Cs2{d, a)Csi{d, a) Pk + Asi{d, a)R + Bsi{d, a)\k 



{Sf[g.g]+n 
Here we set (•S'g^'' [0, f , v] g 



S'i[9.9] + Tlk)^,^, = V2 



Pk- 



R 
1 



Sfi[g,g] + nk) _ = V2 



Pk-2^ 

d , 



(4.16) 



Now we are in a position to write down the RHS of the flow equation with g = g. We 
shall denote it Sk{R) in the following. In Sk{R) we need the inverse operators [T^^^ +71^)^^ 
and {S^£ +nky\ The inversion is carried out in appendix P 1| . Inserting the inverse 
operators into Sk{R) leads to 



SkiR) 



Tr 



(25T2) 



[Pk + AT{d)R - 2\k) Vj + Tr'(iy) \{Pk + Ay{d, a)R - 2a\k) V 



(Pk + AsMR - 2Afe) ' {Pk + Asi{d, a)R - 2a\k) 



-1 



X • 



[FsM a)Pk + As,{d, a)R - 2{a + l)Xk) U 



+Fs2{d,a)JPkJPk 



R 
d-l 
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-2Tr 



(IT) 



E 

Nk i=o 



2Z 



Here we set 



Di{d,0) 



2Tr' 



(0) 



Pk 



R 



dt\z^kk^R^''\Ai{d,0)/k^ 



Ai{d, 0) + km(o\Ai{d, 0)/P) + Asi{d, a)R + BsM 



• (4.17) 



1 - 2^^(^) 



where 



y]N{k) 



-df In Z 



Affc 



(4.18) 



(4.19) 



is the anomalous dimension of the operator / d'^Xy/gR and the prime at denotes the 
derivative with respect to the argument. Furthermore, the new A's and F's introduced 
above are again functions of d and a, tabulated in appendix 0. 

In eq. (|4.17| ) we refined our notation concerning the primes at the traces. From now 
on one prime indicates the subtraction of the contribution from the lowest eigenvalue, 
while two primes indicate that the modes corresponding to the lowest two eigenvalues 
have to be excluded. 

The next step is to extract the contributions proportional to / d'^x^ and / d'^x^/gR 
by expanding Sk{R) with respect to R or r, respectively. Since / d'^x^Jg oc r"', / d'^x^JgR oc 



"'^ ^, only terms of order r'^ and r'^ ^ are needed. This leads to 

1 



5fc(i?) = Tr(2ST2) 



{Pk - 2A 



+ Tr; 



(IT) 



Tr" 



P, - 2a\i ' 



Af 
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-2Tr 



(IT) 



2%) [Pk'Afo 



-i?|AT(rf)Tr(25T2) \{Pk-2X 



+Asi{d,a)Ti'l^) \(Pk - 2aAfc)~V 



(0) 



2t 



<d-2\ 



(4.20) 



Here 0{r'^'^~^) means that terms oc r" with powers n < d — 2 are neglected. 

The term in eq. ( |4.2(J| ) proportional to 6d^2 arises from the last term in eq. ( [4.171 ). 
Contrary to the other terms of eq. ( 4.1?1 ), its expansion does not contain rf-dependent 



powers of r, but is of the form J2 &2m^ with {62m} a set of r-independent coefficients. 

m=0 

As for comparing powers of r, this has the following consequence. Since, for all m > and 
d > 0, —2m = d — 2 is satisfied only if m = and d = 2, and since —2m = d cannot be 
satisfied at all, this term contributes to the evolution equation only in the two-dimensional 
case. Using eq. ( |4.7] ) the piece contributing, i.e. 62m=o may be expressed in terms of 
the operator / d'^x^/gR which yields the last term in eq. (|4.2CI|) . 

The traces appearing in eq. ( [4.20| ) are evaluated in appendix p 2| using heat kernel 
techniques. Then combining the result with the LHS of the fiow equation, eq. ( [4.5| ), and 
comparing the coefficients of the invariants / d'^x^Jg and / d'^x^R leads to the desired 
system of coupled differential equations for Z^k and A^. We obtain 

dt {ZNkXk) = {AK'Y\A7i)-''/Vl^d{d - 1) $i/2(-2Afc/A;2) + d<!>l/,{-2a\k/k^) (4.21) 

-2A,/A;2) + rf$^/2(-2aA,./A;2)l - 2d$i/2(0) 



-d{d 



l)*d/2( 



dtZ^k = -(2/t2)-i(47r)-'^/2A;'^-2|ci(rf) ^\/,_,{-2Xk/e) + C2{d) ^\,^^^{~2a\k/k' 



+C3{d) ^l/^{-2X,/e) + c^{d, a) ^l^^{-2aXk/k^ 



-2^N{k) 



ci{d) $i/2_i(-2A,./A;2) + C2{d) ^\,^^^{-2a\k/k 
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d,2 



2Afc/P 1 - 2a\klk'^ 



(4.22) 



Here $5^^, $^ are cutoff-dependent "threshold" functions defined as 



W 



hi?(0)(y)- 



n > 




n 



r(n) J (y+R{o)(y)+«,)'' 



n > 




(4.23) 



The coefficients q are given by 

- 2£ - lid - 12 



ci{d) 
Ci{d^ a) 



i2(rf- 1; 

ad{d — 2) 



C2{d) 



d^-6 



6d 



, C3(c?) 



rf3 - 4rf2 + 7d - 8 
2(rf- 1) 



rf- 1 



C5(ci) 



2(rf+ 1) 



(4.24) 



d ' ' ' d 

In eq. ( [4.22| ) the terms proportional to 5rf^2 arise not only from the last term of eq. 
( [4.2U| ), but also by evaluating the "primed" traces, i.e. by subtracting the contributions 
coming from unphysical modes, see appendix |B 2| for details. All these contributions are 
obtained by expanding various functions f{R) with respect to R and retaining only the 
term of zeroth order, /(O). As we argued above, these are the only pieces of / which 
may contribute to the evolution in the truncated parameter space. Furthermore, the heat 
kernel expansions of the traces corresponding to differentially constrained fields introduce 
additional contributions oc 6d^2 into eq. ( [4.221 ). 

In appendix P we concentrate on the 4-dimensional case and compare our result for 



Sk{R) and for the corresponding RG flow of Z^k and with the one of ref. ||2T[ where a 
cutoff of type B is used, too. 
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D. The system of flow equations for and 



Now we introduce the dimensionless, renormalized Newton constant 



gk = k''-^Gk = k''-'Z^iG 



d-2 ry-\ 



and the dimensionless, renormahzed cosmological constant 



(4.25) 



(4.26) 



where denotes the corresponding dimensionful, renormahzed Newton constant at scale 
k. Inserting eq. ( |4.26| ) into dt^Z^k^k) leads to the relation 



dtXk = - {2 -r]N{k))Xk + 3271 gkKH-''dt{ZNkXk) 



(4.27) 



Then, by using eq. ( |4.21| ), we obtain the following differential equation for the dimen- 
sionless cosmological constant: 



dtXk = (3x{Xk, gk] a, d) = Ai{\k, gk] a, d) + r]N{k) A2{\k, gk] a, d) 



(4.28) 



The /3-function f3x contains the quantities Ai and A2 which are defined as 
A^{Xk,gk;a,d) = -2\k + {4^)'-"^^ gk^d{d - 1) ^l/^{-2\k) 



+2(i$i/2(-2«Afe)-4(i$i/2(0) 
A^iXk, gk, a, d) = Xk- (47r)i-'^/2 _ 1) ^1^^^{-2X,) + d^\/^{-2aXk) \ . (4.29) 



The corresponding /3-function for gk may be determined as follows. Taking the scale 
derivative of eq. ( |4.25| ) leads to 



dtgk = (3JXk, gk, a,d) = [d - 2 + r]N{k)] g^ 



(4.30) 



For the anomalous dimension riN{k) we obtain from eq. (|4.22| ) 



VN{k) = gk Bi{Xk] a, d) + 7]N{k) gk B2{Xk] a, d) 



(4.31) 
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with Bi, B2 functions oi Xk, d and a given by 



d/2- 



1 1 



+3(5, 



(i,2 



1 - 2Ai. 1 - 2aAi 



S2(Afc; «, rf) = -2(47r)i-'^/2|^i(c/) $;^/2-i(-2Afc) + C2(rf) $i/2-i(-2aAA 
+c^{d) $^/2(-2Afe) + C4(ci, a) $^/2(-2aAfc) 



(4.32) 



.l-2Afc l-2aAfc_ 

Eq. ([4.31|) may now be solved for the anomalous dimension in terms of A^, Qk-, ot and d: 

Qk Bi{Xk;a,d) 



(4.33) 



1 - 9kB2iXk;a,d) ' 

The system of coupled flow equations ( ^4.28| ), ( ^4.30| ) is the main result of this section. 



E. Comparing the cutoffs A and B 



Let us compare the flow equations ( ^4.28| ), ( [4.30| ) of the present paper with those 
obtained in ref. P|. Ref. P] covers the case a = 1 only, and the cutoff used there has a 
different structure than the present one. In ref. 0, the A^S for the cutoff A is formulated 
at the level of the complete field h^^, i.e., symbolically, A^S oc / h^yTZkh^" , while the 
cutoff B of the present paper is based upon a similar action for the component fields: 
A^S oc / hlJRkh^^^^ + / e;.7^fce + ■ ■ ■■ 

For a = 1 the new /3-function for A^, (3^^, agrees perfectly with the result in [0, whereas 
the coefficients B2 in the /3-function for g^, f3g, do not coincide with the corresponding 
results derived there. However, with both cutoffs these coefficients are of the form 

B,{X,; a = l,d) = i(47r)i-'^/2fc'^-2|ei(rf) <l>i/2-i(-2Afc) + e^id) $^/2(-2A,) 
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B2iXk;a = l,d) = -^(47r)i-'^/2A;'^-2|gi(f/)$^/2_i(-2A,) + e2(rf)$^/2(-2Afc)| . (4.34) 

In the present paper (cutoff B) the coefficients are obtained as 

, „ d^- 13rf2 - 2Ad +12 , - 2d^ ~(P -Ad + 2 

= — w^) — ' "'^ ^ = — w^) — ' 

(i^ — 6 , „ d+1 , ^ 

e3(ci) = -4^— , e4(ci) = -24^- (4.35) 

while in (cutoff A) they are given by 

ei{d) = d{d + 1) , e2{d) = -6d{d - 1) , e3{d) = -Ad , ei{d) = -24 . (4.36) 



Upon subtracting the coefficients in eq. ( [4.35[ ) from those in eq. ( |4.36| ) we obtain Aei 



— Ae2 = 12{d'^ + 2d — l)/{d{d—l)) and Aes = — Ae4 = —24/d. Quite remarkably, the sum 
of the deviations Acj vanishes not only in total but also separately for the gravitational 
contributions, involving Ci and 62, and the contributions from the ghost, which contain 
63 and 64. Obviously, this amounts to a shift from the gravitational {p = 1 , n = d /2 — 1)- 
sector to the {p = 2,n = d/2)-sectoT as well as to a shift between the corresponding 
ghost-sectors. The simplicity of this result is somewhat mysterious. 

V. THE FIXED POINTS 

A. Fixed points, critical exponents, and nonperturbative renormalizability 

Because of its complexity it is impossible to solve the system of flow equations for 
and Afc, eqs. ( [4.2(j| ) and ( [4.30D , exactly. Even a numerical solution would be a formidable 
task. However, it is possible to gain important information about the general structure 
of the RG flow by looking at its fixed point structure. 
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Given a set of /3-functions corresponding to an arbitrary set of dimensionless essential 
couplings gi{k), it is often possible to predict the scale dependence of the couplings for very 
small and/or very large scales k by investigating their fixed points. The fixed points are 
those points in the space spanned by the gj where all /3-functions vanish. (The essential 
couplings are those combinations of the couplings appearing in the action functional that 
are invariant under point transformations of the fields.) Fixed points are characterized 
by their stability properties. A given eigendirection of the linearized flow is said to be 
UV or IR attractive (or stable) if, for k —* oo or k 0, respectively, the trajectories are 
attracted towards the fixed point along this direction. The UV critical hypersurface in 
the space of all couplings is defined to consist of all trajectories that run into the fixed 
point for k —>■ oo. 

In quantum field theory, fixed points play an important role in the modern approach 
to renormalization theory P]. At a UV fixed point the infinite cutoff limit can be taken 



in a controlled way. As for gravity, Weinberg |T6| argued that a theory described by a 
trajectory lying on a finite- dimensional UV critical hypersurface of some fixed point is 
presumably free from unphysical singularities. It is predictive since it depends only on a 
finite number of free (essential) parameters. In Weinberg's words, such a theory is asymp- 
totically safe. Asymptotic safety has to be regarded as a generalized, nonperturbative 
version of renormalizability. It covers the class of perturbatively renormalizable theories, 
whose infinite cutoff limit is taken at the Gaussian fixed point g*j = 0, as well as those 
perturbatively nonrenormalizable theories which are described by a RG trajectory on a 
finite-dimensional UV critical hypersurface of a non-Gaussian fixed point g^^j 7^ and are 



nonperturbatively renormalizable therefore K 



Let us now consider the system of differential equations 

k^kg^{k)=(3,ig) (5.1) 

for a set of dimensionless essential couplings g(fc) = {gi(fc), . . . ,g„(fc)}. We assume that 
g,, is a fixed point of eq. (|5.1| ), i.e. /3j(g*) = for alH = 1, . . . , n. We linearize the RG 
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flow about g^, which leads to 

n 

kdkgiik) =J2Bij {gj{k) -g^j) (5.2) 



1=1 



where Bij = c}j/3j(g*) are the entries of the stability matrix B = [Bij). Diagonalizing 

B according to S'^-'^BS' = — diag(6'i, . . . , 6'„), 5* = {V^, . . . , V^), where is the right- 
eigenvector of B with eigenvalue —6i we have 

n 

J2B,,V/ = -9,V/ , J = l,...,n. (5.3) 
i=i 

The general solution to eq. ( |5.2| ) may be written as 

n /k \ 

&(fc)=g« + E^/^/ X • (5-4) 



7=1 \ ^ 



Here 



C/ = E(^"%g.(M (5.5) 
i=i 

are arbitrary real parameters and /cq is a reference scale. 

Obviously the fixed point (7* is UV attractive (i.e. attractive for k 00) only if all Cj 
corresponding to negative 61 < are set to zero. Therefore the dimensionality of the UV 
critical hypersurface equals the number of positive 6j > 0. Conversely, setting to zero all 
Cj corresponding to positive Oj, g^, becomes an IR attractive fixed point (approached in 
the limit k 0) with an IR critical hypersurface whose dimensionality equals the number 
of negative 6*/. 

In a slight abuse of language we shall refer to the ^/'s as the critical exponents. 

Strictly speaking, the solution (|5.4| ) and its above interpretation is valid only in such 
cases where all eigenvalues —6j are real, which is not guaranteed since the matrix B is not 
symmetric in general. If complex eigenvalues occur one has to consider complex C/'s and 
to take the real part of eq. ( ^.4|) , see below. Then the real parts of the critical exponents 
determine which directions in coupling constant space are attractive or repulsive. 



43 



At this point it is necessary to discuss the impact a change of the cutoff scheme has 
on the scahng behavior. Since the path integral for depends on the cutoff scheme, 
i.e. on the A^S chosen, it is clear that the couplings and their fixed point values are 
scheme dependent. Hence a variation of the cutoff scheme, i.e. of TZk, induces a change 
in the corresponding B-matrix. So one might naively expect that also its eigenvalues, 
the critical exponents, are scheme dependent. In fact, this is not the case. According 
to the general theory of critical phenomena and a recent reanalysis in the framework of 
the exact RG equations any variation of the cutoff scheme can be generated by a 
specific coordinate transformation in the space of couplings with the cutoff held fixed. 
Such transformations leave the eigenvalues of the B-matrix invariant, so that the critical 
behavior near the corresponding fixed point is universal. The positions of fixed points are 
scheme dependent but their (non)existence and the qualitative structure of the RG flow 
are universal features. Therefore a truncation can be considered reliable only if it predicts 
the same fixed point structure for all admissible choices of TZk- 

In the context of the Einstein-Hilbert truncation the space of couplings is parametrized 
by gi = A and g2 = g. The /3-functions occuring in the two flow equations 

dAk = f^xi^k, 9k) , dtQk = /3„(Afc, Qk) (5.6) 



are given in eqs. ( [4.28| ) and ( [4. 301) , respectively. As we shall see in subsection B, they 



have a trivial zero at X^, = g^, = 0, referred to as the Gaussian fixed point. The analysis 
of subsection C reveals that there exists also a non-Gaussian fixed point at A* ^ 0, 
g^. 7^ 0. In subsection C we study its cutoff dependence and the cutoff dependence of the 
associated critical exponents employing the above /3-functions of type B as well as those 
of refs. [p2| , p!9| based on the cutoff type A, with the families of shape functions 



or 



^ inserted. 
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B. The Gaussian fixed point 



In this subsection we discuss the features of the Gaussian fixed point (A*, g^) = (0, 0). 
In order to investigate the RG flow in its vicinity we expand the /3-functions in powers 
of A and g according to eq. (H5) of appendix ^ and read off the B-matrix. It takes the 
form 

B = ' . (5.7) 

\^ d~2 j 

Here z/^ is a d-dependent parameter defined as 

j.,^(rf_3)(4^)i-|$i^^(0). (5.8) 

Diagonalizing the matrix (|5.7| ) yields the (obviously universal) critical exponents 6i = 2 
and 02 = 2 — d which are associated with the eigenvectors = (1,0)'^ and V'^ = 
{vdi 1)^- Hence, for the linearized system obtained from ( |H5| ) the solution ( [5 .41 ) assumes 
the following form: 

M = {>^ko - t^dgko) j ^^d3ko\jrj 

9^ = 9kA-^] . (5.9) 



Since the expanded /3-function (3g of eq. ( P5| ) is Afc-independent up to terms of third 
order in the couplings we can easily calculate also the next-to-leading approximation for 
gk near the fixed point. In terms of the dimensionful quantity this improved solution 
reads 

-1 



Gk — Gk. 



l-ujaGkAkt'-k''-'] (5.10) 



with 



1 ^ /n ,^ o/. M-^ f(^ + 2)(c?^-6rf2 + 3d-6) ^1 

- ( d[d-m-2) + H ^%M\ (5.11) 
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a d- and a-dependent parameter. For k <^ Icj^G^pl^^/*^'^^^) and with the reference scale 
kf) = (which is admissible only for trajectories lying on the IR critical hypersurface of 
the fixed point) eq. ( |5.10| ) yields 

Gk = Go[l- L^dGok'"-' + O {Gie^^-^^)] . (5.12) 

For the dimensionful cosmological constant we obtain from eq. (|5.9| ) 

Xk = Xko + i^dGko{k''-k^) . (5.13) 

Apart from the different expression of Ud due to the new cutoff, the solutions ( |5.12| ) 
and ( p.l3| ) coincide in = 4 dimensions with those derived in ref. by using a similar 
approximation scheme, see eqs. (5.18) and (5.25) of this reference. 

Let us now analyze the scaling behavior near (A*,^'*) = (0,0). Since 6'i = 2 > 
the V^^-eigendirection, which coincides with the A-direction, is IR repulsive (and thus UV 
attractive). For d < 2, $2 is positive which implies that the Gaussian fixed point is UV 
attractive for any direction in the 2-dimensional parameter space. 

For d > 2 we have < so that the ^^-eigendirection is IR attractive (and UV 
repulsive). Hence, in this case both the UV and the IR critical hypersurface of the 
Gaussian fixed point are one- dimensional, i.e. they consist of a single trajectory. For the 
IR critical trajectory that hits the fixed point in the limit — we have 

Xk = UdGkk'^ Xk = i^d9k (5.14) 

for sufficiently small values of k with Gk given by eq. ( |5.12| ). Since $^/2(0) depends 
on the shape function R^^\ Ud is not a universal quantity. Therefore the slope of the 



distinguished trajectory ( 5.14 ) is not fixed in a universal manner. This is in accordance 
with the general expectation that the eigenvalues of B should be universal, but not its 
eigenvectors. 

For d 2 the parameters ^d/2~ii^) ^"^^ ^d/2i^) appearing in Ud are scheme dependent 
as well. Furthermore, Ud is a function of the gauge parameter a. Hence Ud is a nonuniversal 
quantity, too. In the most interesting case of d = 4 dimensions it takes the form 
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UJ4 



24:71 



13$J(0) + (55 + 24a)$2(0) 



(5.15) 



Since $i(0) and $2(0) positive for any admissible shape function we can infer from 
eq. ( ^.15|) that is positive for all a > ao and negative for all a < aQ. Here = 
[— 13$]^(0)/<l'2(0) ~ 55]/24 is a negative number of order unity. Thus, if we identify 
Einstein gravity with the theory described by the IR critical trajectory of the Gaussian 
fixed point, eq. ( |5.12| ) implies that Einstein gravity is antiscreening for all a > aQ, i.e. Gk 
decreases as k increases. On the other hand if a < ao gravity would exhibit a screening 
behavior. As argued in subsection [IV A| , the gauge parameter should be regarded as a 
scale dependent parameter in an exact treatment where it is expected to approach the 
fixed point value a* = 0. Setting a = from the outset we may conclude that the physical 
Gk displays the antiscreening behavior found for a > ao- 

In ref. a similar result was obtained with a cutoff of type A, while the above 
calculation employed the cutoff B. The only difference between our result for the behavior 



of Gk and the one obtained in [g^ lies in the slightly differing value of ao which is a 
scheme dependent parameter. For the cutoff type A, 004 = 00^^^ is given by (see 



(A) 



1 



;i8 + 6a)$^(0) - ^^(O) 



(5.16) 



so that ao = [$j;(0)/$2(0) — 18]/6 in this case, while z/4 is the same with both cutoffs. For 
a = 1, eq. boils down to lu^^^ = [24<l>^(0) - $^(0)] /(Gtt) which equals the result 

obtained in the original paper [||. This is because for a = 1 the cutoff type A coincides 
with the one used in ref. [^]. For a comparison of this result with the one for the cutoff 
type B we insert a = 1 into eq. ( p7[5| ) which yields cjf ^ = [13$}(0) + 79$2(0)] /(247r). 
Using the exponential shape function R^'^^ with s = 1 we have $}(0) = vr^/6, $2(0) — 1 



so that 00^^ ^ 1.33 for the cutoff type B, which lies rather close to the value uj^ 



(A) 



1.19 



obtained in P,|22[] for the cutoff type A. Furthermore, we have $2(0) — 2C(3) where ( 
denotes the zeta function, and thus ^ 0.19 with both cutoffs. 
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C. The non-Gaussian fixed point 



Now we turn to the nontrivial zeros of the set of /3-functions {f3x,f3g} given by eqs. 
( [4.28| ), ( [4.301) . Such non-Gaussian fixed points {\*,g*) 7^ (0,0) satisfy the condition 



r]N, = 2-d (5.17) 
which follows immediately from eq. ( [4.30|) . 



1. In 2 + e dimensions 

As a warm up we consider the case of d = 2 + e dimensions with < |£:| ^ 1 which 
can be dealt with analytically and for which the existence of a non-Gaussian fixed point 
has already been shown p6|j6|,^^-|3^ ■ this case the condition ( p.l7|) takes the form 



Vn* = —£ with 



l-^7*(^)i?2(A.(e);a,2 + 5) " ^"'''^ 



Solving eq. ( p. IS] ) for g^{e) and expanding the result with respect to e leads to 

g4e) = -[Bi{X40y,a,2)r' e + 0{e^) . (5.19) 

Furthermore, expanding also Px{X^,{e), g^,{e); a, 2 + e) with respect to e and equating equal 
powers of e yields 

K{e) = [Si(0; a, 2)]"' $^(0) e + O (e^) . (5.20) 

In particular we obtain A*(0) = which implies g^:{e) = — [Bi{0; a, 2)]~^ e + O [e"^). 

The parameters -Bi(0; a,2), i = 1, 2, may be obtained from the zeroth order terms of 
the expansions Bi{Xk; a,2 + e) = Bf\Xk; a) + Bl_^\Xk; a)e + 0{6^) which take the form 

B['\Xk; a) = -f (1 - 2Afe)-^ - ^ (1 - 2aA,)-' + 4$?(-2A,) + 6$?(-2aA,) - y , 
5f (Afc; a) = y (1 - 2Afc)-' + ^ (1 - 2«Afc)-' - 2¥,{-2Xk) - 3^{-2aXk) . (5.21) 
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Inserting = into eq. (|5.21|) and using that $i(0) = 1 is scheme independent yields 



5i(0;«,2) = 5fH0;a) = -y , 

^2(0; a, 2) = Bi°\0; a) = 6 - 5$?(0) . (5.22) 

In contrast to the universal quantity Bi'\0;a), B2^\0;a) depends on the shape of R^^^ 
via $i(0). However, i?2°^(0;Q;) does not enter the leading order term of X^{e) and g^{e), 
which may now be written as 

3 . .^/_2 



9.{e) = -e + 0[e'} . (5.23) 

The leading order term of X^{e) is nonuniversal since it contains the scheme dependent 
parameter $i(0). This is not the case for g^{e) whose leading order contribution has a 
universal meaning. 

Let us now analyze the scaling behavior near the non-Gaussian fixed point ( 5.23 ). One 
finds that the associated B-matrix is of the form 
/ 



B 



V 



-2 + £ + O (e^) -2$}(0) + O(e) 

O (e^) -e + O (e^) 



(5.24) 



From (|5.24|) we obtain the critical exponents 9i = 2 — -"^^^ e + O (e^) and 62 = £ + 



Ole"^). 61 and 62 are scheme independent up to terms of 0{e'^). 61 depends on the 
gauge parameter. For a = 1 the critical exponents coincide with those following from 
the /3-functions of ref. 0, where the cutoff type A is used. These findings nicely confirm 
that, to lowest order, the critical exponents are the same for the cutoffs A and B and are 
independent of R^^\ 

For £ > both critical exponents are positive. Hence the non-Gaussian fixed point 
( ^.23| ) is UV attractive for all trajectories so that the condition for the asymptotic safety 
scenario is met. It is interesting to investigate whether this result stabilizes in the sense 
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that more general truncations including higher powers of the curvature tensor reproduce 
this fixed point and lead to a finite-dimensional UV critical hypersurface. In we will 
discuss this point in detail. 



2. Location of the fixed point {d = 4) 

In d = 4 dimensions, and for the cutoff A, the non-Gaussian fixed point of the Einstein- 
Hilbert truncation was first discussed in |17], and in ref. [19| the a- and /^'-''-'-dependence of 



its projection (0,(7*) onto the gf-direction has been investigated. However, since for a^l 
the cutoff of type A is introduced by an ad hoc modification of the standard one-loop 
determinants it is not clear whether it may be expressed in terms of an action A^S", except 
for the case a = 1 0. Since a specification of A^S" is indispensable for the construction 
of Ffc, the status of the results derived in [|TU| is somewhat unclear. In the following we 



determine the fixed point properties using different cutoffs of type B, for which a A^S* is 
known to exist, and compare them to the analogous results for the cutoff A. 

In a first attempt to determine the non-Gaussian fixed point we neglect the cosmolog- 
ical constant and set = A* = 0, thereby projecting the renormalization group fiow onto 
the one-dimensional space parametrized by g. In this case the non-Gaussian fixed point 
is obtained as the nontrivial solution of f3g{0,g^:;a,d) = 0. It is determined in appendix 
0with the result given by eq. (p^). In order to get a first impression of the position of 
we insert the exponential shape function with s = 1 into eq. ( |H2|) and set d = 4, a = 1. 
We obtain g^, ^ 0.590. 

Assuming that for the combined X-g system both g^ and A* are of the same order 
of magnitude as g^ above we expand the /3-functions about {Xk,gk) = (0,0) and neglect 
terms of higher orders in the couplings. Again in appendix ^ we determine the non- 
Gaussian fixed point for the corresponding system of differential equations. Inserting the 



shape function (|3.7| ) and setting c? = 4, a = 1, we find {X^,g^) ^ (0.287,0.751). 



In order to determine the exact position of the non-Gaussian fixed point (A*,^'*) we 
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20 20 

(a) (b) 

FIG. 1. (7# as a function of s and a from the approximation A/; = A,, = 0, using (a) the cutoff 
type A, and (b) the cutoff type B, with the family of exponential shape functions (|3.8|) inserted. 



have to resort to numerical methods. Given a starting value for the fixed point, e.g. one 
of the approximate solutions above, the program we use determines a numerical solution 
which is exact up to an arbitrary degree of accuracy. Under the same conditions as above, 
i.e. s = 1, c? = 4, we obtain 

f (0.348,0.272) for a = 1 
iK,g.) = {' . (5.25) 

[ (0.339, 0.344) for a = 

Next we study the gauge and scheme dependence of the non-Gaussian fixed point. 
The scheme dependence is investigated by looking at the s-dependence introduced via the 
family of exponential shape functions ( |3.8D where s parametrizes the profile of R^^\ 

FIG. |l] shows g^{a,s) obtained from the approximation Xk = K = 0, while FIGS. |] 
and D display the (exact) functions g^{a, s) and A* (a, s) resulting from the combined \-g 
system. In each of these figures the plot on the LHS (i.e. FIGS. 1,2, 3(a)) is obtained 
from the cutoff type A and the one on the RHS (i.e. FIGS. 1,2, 3(b)) is obtained from the 
cutoff type B used in the present paper. 

Our results establish the existence of the non-Gaussian fixed point in a wide range 
of a- and s-values. As expected, the position of the fixed point turns out to be s-, i.e. 
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20 20 

(a) (b) 

FIG. 2. The exact (7* as a function of s and a from the combined X-g system, using (a) the 
cutoff type A, and (b) the cutoff type B, with the family of exponential shape functions ( |3.8| ) 
inserted. 

scheme dependent, but the crucial point is that it exists for any of the cutoffs employed. 
This is one of the important results of our analysis because it gives a first hint at the 
reliability of the Einstein-Hilbert truncation. 

As for the a-dependence, a = is, in principle, the only relevant case since according 
to subsection [IV A| , a = is assumed to be the physical value of the gauge parameter. In 
practical calculations a = 1 is often used instead, because this simplifies the evaluation of 



the flow equation considerably. In [Q, for instance, all calculations are performed with 
a = 1 for this reason. Therefore it is necessary to compare the gauges a = and a = 1 in 
order to judge whether the results obtained by using a = 1 are a sensible approximation 
to the physical case a = 0. Here we see that this is indeed the case. 

As for comparing different types of cutoffs, we recognize from FIG. ^ that, in the 
approximation = A* = 0, the s- dependence of g^, is much weaker for type B than 
for type A. Contrary to this, both cutoffs yield nearly the same results for (7* and A* 
if we consider the combined X-g system, see FIGS. ^ and |^. Furthermore, the scheme 
dependence of g^, in FIG. ^ is stronger than in FIG. |I|(b), but much weaker than in FIG. 



|l](a). FIG. |T](a) reproduces the result of ref. |T9| obtained from the cutoff A, see FIG. 2 
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(a) (b) 

FIG. 3. The exact A=k as a function of s and a from the combined A-i^ system, using (a) the 
cutoff type A, and (b) the cutoff type B, with the family of exponential shape functions (|3.8| ) 
inserted. 

of this reference. 

It should be noted that we are forced to restrict our considerations to shape functions 
( p.8|) with s > 1. This is because for s < 1 the numerical integrations are plagued by 



convergence problems which is due to the fact that in d = 4 dimensions the threshold 
functions in fB^^ and (3g diverge in the limit s — > 0, see also |T^. 

Because the scale k enters the flow equation via TZ^ as a purely mathematical device 
it is clear that the functions k ^ Xk, gk and their UV limits A,,, g^, are scheme dependent 
and not directly observable therefore. It can be argued that the product (^^A* must be 
scheme independent, however. While k and, at a fixed value of k, Gk and A^ cannot be 
measured separately, we may invert the function k \—>- Gk and insert the result k = k{G) 
into Afc. This leads to a relationship between Newton's constant and the cosmological 
constant which, at least in principle, could be tested experimentally: A = X{G). In 
general this relation depends on the RG trajectory chosen (specified by its IR values Aq 
and Go, for instance), but in the fixed point regime all trajectories approach A^ = A*/c^ 
and Gk = g*/k'^ which gives rise to 

X{G) = ^ . (5.26) 

Gr 

Eq. (|5.26|) is valid if A ^ rripy and G <^ nipy. (We define the Planck mass in terms of 
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— 1/2 ~ 

the IR limit of Gk, ra-pi = Gq .) Assuming that A and G have the status of observable 
quantities, eq. (|5.26| ) shows that (^^.A* must be observable, and hence scheme independent, 
too. (For a related discussion see PB|.) Below the Planck regime the function A(G) 
becomes much more complicated than A oc 1/G (which follows already from dimensional 
analysis) because the dimensionful quantities Aq and Gq enter explicitly there. 



0.4 



0.3 



0.2 



0.1 




FIG. 4. (a) s-parametric plot of (A=„(s), (7*(s)) in the range 1 < s < 50 for various values of 
a. Each curve starts on the left at s = 1, and ends on the right at s = 50. (b) a-parametric 
plot of (A*(a), (7*(a)) for various values of s. In both (a) and (b) the cutoff type B is used, with 
the family of exponential shape functions ( |3.8D inserted. 



As for the universality of (?*A*, it is also interesting to note that, for any k, the product 
gk^k = GkXk is essentially the inverse of the on shell value of F^. The stationary points of 
( |1.2| ) with (jf^i, = g^y satisfy Einstein's equation G^y = —XkQ^iv Hence R = 4Afc, so that 



from ( |1.2| ) in four dimensions, 

Ffc[on shell] 



(5.27) 



SirGkXk S-KgkXk 

Here we used that, for dimensional reasons, / d^x^/g = v/\\ where f is a finite, positive 
constant for any solution with a finite four-volume. 

Quite remarkably, the universality of the product gf^A* is confirmed by our results in 
a rather impressive manner, as is illustrated in FIGS. FIG. ^(a) contains several 
parametric plots of (A*(s), (7*(s)) for various values of a, obtained from the /3-functions 
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FIG. 5. g^, X^, and g^=A* as functions of s for (a) a = 0, 1 < s < 30, (b) a = 0, 1 < s < 5, 
(c) a = 1, 1 < s < 30, and (d) a = 1, 1 < s < 5, using the cutoff type B with the family of 
exponential shape functions (|3.8|) inserted. 



( [4.28| ) and (|4.30|) which are based on the cutoff type B. The hyperbolic shape of these 
plots is a first hint at the s-independence of the product (^^.A*. Its direct confirmation 
is supplied by FIG. |] which shows g^, A*, and (^^.A* as functions of s for a = (FIG. 
^(a),(b)) and a = 1 (FIG. ||(c),(d)), again using the cutoff type B. In FIG. |^(a),(c) these 
functions are plotted in the range of values 1 < s < 30 while FIG. ^(b),(d) contains the 
sector corresponding to 1 < s < 5 where the largest changes in A^, and occur. In any of 
these figures the product of A* and is almost constant for the whole range of s-values 
considered. Its universal value is 
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g» A.,: 

1.5 2 2.5 3 3.5 4 4.5 5^^ 

(b) 

FIG. 6. 5*, A*, and g^X^: as functions of s for (a) a = 0, and (b) a = 1, using the cutoff type 
A with the family of exponential shape functions (^^) inserted. 



(7* A* 



0.12 for a = 1 
0.14 for a = 



(5.2^ 



Obviously the difference between the physical case a = and the case preferred for 
technical reasons, a = 1, is rather small. 

It is reassuring to see that employing the /3-functions of refs. [22,1^ which are based 
on the cutoff A we obtain almost identical results. They are illustrated by means of FIG. 
H which shows g^,, X^,, and (yf^A* as functions of s, 1 < s < 5, for a = (FIG. H(a)) and 
a = 1 (FIG. 1(b)). 

It is also interesting to compare the above results with those obtained from a different 
shape function. FIG. ^ displays the behavior of (7*, A*, and their product resulting from 
the /3-functions ( [4.28[ ) and ( [4.30| ) of the cutoff B, with the family of shape functions 
R^'^\y;b) with compact support, eq. ( |3.9|) , inserted. Here b G [0,1.5) parametrizes the 
profile of these shape functions. In FIG. |^(a) we present a parametric plot of {X^:{b), g^{b)) 
for a = starting at 6 = and ending at 6 = 1.5. Furthermore, FIG. ^(b) shows the 
product (7=kA=k as a function of b for a = and a = 1. For b < 1.2, the parametric plot in 
FIG. 0(a) exhibits an approximately linear behavior which leads to a (7*A*-plateau in FIG. 
|^(b) where (7* A* is nearly constant. Remarkably, the position of these plateaus coincides 
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(a) (b) 

FIG. 7. (a) 6-parametric plot of (A*(6), for a = 0, and (b) g*A* as a function of b for 

a = and a = 1, using the cutoff type B with the family of shape functions with compact 
support ( ^^ ) inserted. 

quite precisely with those of the corresponding plateaus in FIGS. ^ and ^ obtained with 
the other cutoffs. As for the quality of the Einstein- Hilbert truncation this result is rather 
encouraging. 

For b > 1.2 the curves in FIG. ^ have a rather strong and erratic 6-dependence. This is 
because R^^\y;b) approaches a sharp cutoff as b 1.5, which introduces discontinuities 
into the integrands of the threshold functions $^ and Already for b ^ 1.2 the 
/3-functions start to "feel" the sharp cutoff limit so that the results cannot be trusted 
beyond this point. 

As for the a-dependence of the fixed point, in FIG. ^(b) we present parametric plots of 
(A*(a), for various fixed values of s. Here we used the cutoff type B with the shape 

function ( p.8| ) inserted. These plots start at positions in the X-g plane which correspond 
to a = and which are different for the distinct s-values. As a ^ oo all curves run into 
the Gaussian fixed point. 

To summarize: 

a) In 4 dimensions, the Einstein-Hilbert truncation leads to a non-Gaussian fixed point 
with positive values of A* and for all admissible cutoffs, both of type A and type B. The 
scheme independence of this prediction is a nontrivial result. (To emphasize this point we 
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mention that in higher dimensions, where the Einstein-Hilbert truncation is less rehable, 
the fixed point exists or does not exist depending on the cutoff chosen 

b) Universal quantities are strictly cutoff independent only in an exact treatment. Any 
truncation leads to a scheme dependence of these quantities. The extent of this scheme 
dependence is a measure for the reliability of the truncation. The product g^\^ is an 
example of a universal quantity. While we find a considerable scheme dependence of 
g,, and A^, separately, their product is scheme independent at a quite amazing level of 
accuracy, see FIG. |^. As for the reliability of the Einstein-Hilbert truncation, we consider 
this result a highly nontrivial confirmation of our assumption that the region of parameter 
space where the fixed point occurs is well described by this truncation ansatz so that the 
fixed point also exists in the exact theory and is not a truncation artifact. 



3. Higher and lower dimensions 

Before continuing our analysis of the d = 4 dimensional case we study the d-dependence 
of the non-Gaussian fixed point. This is done by means of the parametric plots in FIG. 
H which are obtained from the /3-functions of type B, eqs. (|4.28| ) and ( [4.30| ), with the 
shape function (|3l^ ) with s = 1 inserted. FIG. H(b) shows {X^{d), g^{d)) in 2 < d < 4 for 
a = 1. Remarkably, this plot is almost identical with that in FIG. 4 of ref. which was 
derived from the /3-functions based on the cutoff type A. 

In order to gain information about the behavior of the fixed point in (i > 4 we also 
plotted {X^{d), g^{d)) in 2 < (i < 6, but this time for a = 0, see FIG. ^(a). For d beyond 
d ~ 5, the value of g-t increases significantly, whereas A* seems to approach a constant 
value. 

However, this result might be a truncation artifact since it is plausible to assume 
that in higher dimensions the Einstein-Hilbert truncation becomes less reliable. This 
is because with increasing dimensionality d the number of terms which are relevant at 
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FIG. 8. d-parametricplotof (A*(ci),c/*(d)) for (a) a = 0, 2 < d < 6, and (b) a = 1, 2 < d < 4, 
using the cutoff type B with the exponential shape function ( ^ ) with s = 1 inserted. 



the non-Gaussian fixed point and which are neglected in the Einstein-Hilbert truncation 
increases most probably. At the non-Gaussian fixed point the scaling dimensions of local 
invariants such as / / y^-R, / ^/gR^■, J ^/gR^uR^'^ , etc. are not known a priori. We 
only know that with respect to the Gaussian fixed point all local monomials R, R^, ■ ■ ■ are 
relevant or marginal if their canonical mass dimension does not exceed d. A sensible 
truncation should retain at least the relevant terms, whence it is clear that the number 
of terms needed increases with the dimensionality d. (In 4 dimensions, / y/g and / ^/gR 
are relevant and the (curvature) ^-invariants are marginal.) By analogy we expect that 
the description of the non-Gaussian fixed point, too, requires increasingly high powers of 
the curvature when d is increased [BT . 



4- The critical exponents (d=4) 

Let us now return to the 4-dimensional case and analyze the critical behavior near the 
non-Gaussian fixed point. In order to get a first impression of its features we restrict our 
considerations to the cutoff type B with the exponential shape function ( p.8|) with s = 1 
and to the gauge a = 1. In this case we have (A^,,^',,) = (0.348,0.272), see above. The 
corresponding B-matrix assumes the form 
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B 



/ 



V 



-0.187 5.129 
-3.228 -2.907 



(5.29) 





It leads to a pair of complex critical exponents 6i = 6' + iO" and 02 = 61 = d' — iO" . For 
the real quantities 6' and 9" we find 9' = 1.547 and 9" = 3.835. (In general we define 9i 
as the critical exponent with the positive imaginary part so that 9" > 0.) The behavior of 
Afc and near the fixed point is described by the real part of eq. (|5.4|) in this case. Using 
that V"^ = {V^y, and setting = V and Ci = C, the general solution to the linearized 
flow equation may then be written as 

2<j [ReC cos(^"t) + ImC sin(^"t)]Rel/ 

+ [ReC sin(^"t) -ImC cos (^"t)] Iml^je"^'* . (5.30) 

Here t = ln(fc/fco)- Obviously the non-Gaussian fixed point is UV (IR) attractive if 
9' = Re^i = Re^^2 > (< 0). The imaginary parts ±9" of the critical exponents do not 
infiuence the stability of the fixed point. They only give rise to a rotation of the vector 
(Afc — A*, gk — g^)^ about the fixed point. 

In the case under consideration we have 9' > which implies that the non-Gaussian 
fixed point is UV attractive in both directions of (A, g)-spa.ce. All RG trajectories which 
reach its basin of attraction spiral into the fixed point for k ^ oo. Thus, the Einstein- 
Hilbert truncation predicts all the ingredients which are necessary for the asymptotic 
safety scenario and the nonperturbative renormalizability of 4-dimensional quantum grav- 
ity. Clearly the dimensionality of the UV critical hypersurface cannot be determined 
within the present approach. We shall come back to this question in the framework of a 
more general truncation including higher derivative terms . 



As discussed in subsection |V A| the critical exponents are universal in an exact treat- 



ment, in contrast to g^ and A^,. However, in a truncated parameter space a scheme 
dependence is expected to occur as an artifact of the truncation. Therefore we may use 
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FIG. 9. (a) 9' = Re 01, and (b) 9" = Im^i as functions of s for various values of a, using the 
cutoff type B witli the family of exponential shape functions (|3.8|) inserted. 



this scheme dependence of the critical exponents to judge the quaUty of our truncation. 
Also in this respect the Einstein-Hilbert truncation yields satisfactory results, which we 
display in FIGS. |9|-|TT|. First of all it should be noted that the critical exponents obtained 
from our numerical analysis have a nonzero imaginary part in any of the cases considered. 
FIGS. P|-PT](a) show their real part 6' while their imaginary part 6" is depicted in FIGS. 
|-|ll](b). FIGS. I and ^ are derived from the /3-functions (|]28D and (WM ) based on the 
cutoff B of the present paper, with the family of shape functions ( |3.8[ ) (in FIG. |]) and 
(in FIG. |TU|) inserted. For comparison, FIG. |TT] is obtained from the /3-functions 



of refs. p2| , p!9| derived from the cutoff type A, with the family of shape functions (|3.8|) 
inserted. The figures contain various plots for distinct values of a, which describe the s- 
or 6-dependence of 6' and 6". 

The 6''-plots in FIGS, ^(a) and |ll](a) for the cutoffs B and A, respectively, exhibit a 
similar s-dependence, and the same holds true for the 6'"-plots in FIGS. |^(b) and |ll](b). 
Moreover, in the range 1 < s < 5 one recognizes a total variation of both 6' and 6" 
which is approximately of the order of magnitude of 6' and 6", respectively. For s > 5 
there remains only a rather weak dependence on s, such that the functions 0'{s) and 6"{s) 
develop a plateau-like shape. The s-dependence is slightly stronger for the cutoff type B 
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FIG. 10. (a) 9' = Re^i, and (b) 6" = Im^i as functions of b for a = and a = 1, using the 
cutoff type B with the family of shape functions with compact support ( ^ ) inserted. 

than for type A, and the positions of the "plateaus" are different for both cutoff types. 
For a = 1, for instance, we have 9'{s = 30) ~ 1.56 and 9"{s = 30) ~ 3.06 using type A, 
while employing type B yields 6'{s = 30) ~ 1.75 and 6"{s = 30) ~ 2.31. These differences 
have to be interpreted as truncation artifacts. 

The above results may now be compared with those for the cutoff B with the shape 
functions with compact support, eq. (|3.9|) , inserted. Our results are shown in FIG. 0. 
Remarkably, 6' is almost constant in the range < 6 < 1.2. Furthermore, the change in 
6" is only very weak in this region and occurs mainly for b ^ 0.4. For b ^ 0.4 6" is almost 
constant as well. In the region 1.2 < 6 < 1.5 both 6' and 6" strongly vary with b which 
is caused by the sharp cutoff limit as 6 ^ 1.5, see above. The position of the plateaus 
is nearly identical with those obtained above by using the same cutoff B, but the shape 
function (|3]|). For a = 1 we have e'{b = 0) ^ 1.69 and e"{b = 0) ^ 2.56. 

As discussed above, a = is assumed to be the physical value of the gauge parameter 
while in practical applications a = 1 is often preferred. For a further justification of 
this approximation we compare the critical exponents obtained from a = 1 with those 
obtained from the (physical) gauge a = 0. The results are qualitatively the same, but 
quantitatively we find a relative deviation of about 10 to 20 per-cent. This has to be kept 
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FIG. 11. (a) 9' = Re^i, and (b) 9" = Im^i as functions of s for a = and a = 1, using the 
cutoff type A with the family of exponential shape functions { \i.8\j inserted. 



in mind when calculations are performed with a = 1. 

To summarize: For all admissible cutoffs, both of type A and type B, the non- 
Gaussian fixed point is UV attractive in both directions of parameter space. It is char- 
acterized by a pair of complex conjugate critical exponents which leads to spiral- type 
trajectories in its vicinity. The exact critical exponents are universal. Those obtained 
from the Einstein-Hilbert truncation are approximately scheme independent, but their 
scheme dependence is stronger than that of (/^A*. This might be related to the fact that 
the / y^i?^-term neglected by the Einstein-Hilbert truncation is (very weakly) relevant 
at the non-Gaussian fixed point. We shall come back to this point in ref. p4| . 



VI. THE GRAVITON PROPAGATOR AT LARGE MOMENTA 

The UV fixed point is characterized by an anomalous dimension t] = rj^ig*, A*) = —2. 
We can use this information in order to determine the effective momentum dependence of 
the dressed graviton propagator for momenta in the fixed point regime, i.e. for ^ mpj. 
Expanding the truncated F^ about fiat space and omitting the standard tensor structures 
we find the inverse propagator Qk{p)^^ Zi\i{k)p'^. The conventional dressed propagator 
Q{p) contained in F = F^^q obtains from Qj^ in the limit k ^ 0. Assuming that the actual 
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physical cutoff scale is the momentum itself (which is indeed true for p'^ > ^ ^pi); 
the fc-evolution of Qk{p) stops at the threshold k = . Therefore 

G{j>)-^ oc ZN{k= Vp2^ oc (6.1) 

because Z^lk) oc k~''^ when rj = —dtlnZN is (approximately) constant. In d dimensions, 
and for r] ^ 2 — d, the Fourier transform of G{p) oc l/{p'^Y~^^'^ yields the following 
Euclidean propagator in position space: 

\x — y\ 

This is a standard result well known from the theory of critical phenomena, for instance. 
In the latter case it applies to large distances, while in quantum gravity we are interested 
in the extreme short distance regime governed by the UV fixed point. However, as it 
stands eq. ( |6.2|) is not valid for the case of interest, d = 4 and r] = —2. For rj = —2 
the dressed propagator is Q{p) = which, in c? = 4 dimensions, has the following 
representation in position space: 

g(x; y) = In (/i \x - y\) . (6.3) 

Here fi is an arbitrary constant with the dimension of a mass. Usually in ordinary mat- 
ter field theories on flat space a 1/p^-propagator is considered problematic because it is 
incompatible with a positive spectral density. It is by no means clear, however, that 
the familiar notions of positivity, analyticity, and causality which are appropriate for 
Minkowski space are of relevance to the situation under consideration. For k ^ oo the 
"on shell" spacetimes of Tk have a large curvature oc fc^, after all. 

So let us assume that the improvement ( |6.1| ) is indeed correct and that in the fixed 
point regime, i.e. for distances much smaller than the Planck length, the effective graviton 
propagator has a logarithmic dependence on the distance. This result is quite remarkable 
because it implies a kind of dimensional reduction from 4 to 2 dimensions. In fact, eq. 
( |6.3|) is precisely what one obtains from a standard l/p^-propagator in ci = 2 dimensions. 
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This means that, in a certain sense, spacetime appears to be two-dimensional when it is 
probed by a very high energetic graviton. 

Since, symbohcally, R = ddh, the propagator ( |6.3|) yields the curvature-curvature 
correlation function 

(R(x)R(y))oc-^. (6.4) 

[x-y) 

Its short-distance singularity has to be contrasted with the l/{x — y)^ behavior one finds 
at tree level. Here R stands for the curvature scalar or for any component of the Riemann 
or Ricci tensor. 

Switching for a moment to spacetimes with a Lorentzian signature, it is interesting to 
look at the linearized gravitational field produced very close to a static source. Decom- 
posing X = (x'',x), the relevant Green's function for static problems reads 

oo 

6^stat(x;y)= J dy^g(x\i^;y\y) . (6.5) 

— oo 

In our case this is the 3-dimensional Fourier transform of l/|p|^, i.e. 



^stat(x; y) 



2v72^ 1 



V"V 



y) = -^|x-y| (6.6) 



provided |x — y| ^ nipl. In an, admittedly somewhat naive, Newtonian language this 
result would mean that a point mass located at y = creates a gravitational potential 
which behaves as $(x) oc |x| as long as |x| is much smaller than the Planck length. Prob- 
ably this linear potential is related to a similar phenomenon shown by the renormalization 
group improved Schwarzschild black hole which has been constructed recently ||13|. The 
radial geodesies in this spacetime, in a fully relativistic treatment, experience a lineai|^ 
repulsive "potential" close to the core of the black hole. 



ttin the notation of ref. [^] the potential is linear if 7 = 0. This corresponds to the case 
k (X 1/r which ignores the impact the background curvature has on the cutoff identification, see 
eq. (4.35) of This is consistent with the fact that eq. (|6.3| ) is valid for a flat background. 
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VII. SUMMARY AND CONCLUSION 



On the formal side, the main result of this paper is the construction of a new exact 
RG equation for the gravitational effective average action. It is formulated in terms of 
the component fields appearing in the TT-decomposition of the metric. It is defined on a 
sufficiently large class of background spacetimes so as to facilitate the projection of the RG 
flow onto very general truncated parameter spaces. It also helps in finding admissible IR 
cutoffs. A formalism of this kind is mandatory for truncations including higher-derivative 
invariants, matter fields, or a running gauge fixing, for instance. In a forthcoming paper 
23 we shall use this technology in order to explore a more general truncation with a 



i?^-term. In the present paper we analyzed the Einstein- Hilbert truncation with the new 
equation and a new cutoff. 

After deriving the exact functional equation in section we discussed in section |T| a 
general strategy for constructing an appropriate cutoff operator in the context of a class 
of truncations which is still very general. In section |IV| we specialized for the Einstein- 
Hilbert truncation and derived the nonperturbative ,3-functions which govern the RG 
evolution of and A^. 

On the applied side, most of our results concern the non-Gaussian fixed point of the 
X-g system in 4 dimensions which we analyzed in section |V|. If this fixed point is actually 
present in the exact theory, its importance can hardly be overestimated. Its existence 
would imply that in spite of its notorious perturbative nonrenormalizability quantum 
Einstein gravity is most probably renormalizable at the nonperturbative level and thus 
qualifies as a fundamental (microscopic) quantum theory of gravity. Clearly the crucial 
question is whether the fixed point of the Einstein-Hilbert truncation is indeed genuine 
or merely a truncation artifact. 

In order to get an impression of the reliability of the Einstein-Hilbert truncation we 
investigated how its predictions change when we vary the cutoff which is built into the 
RG equations. We used both the original cutoff of type A, formulated in terms of h^^, and 
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the new cutoff of type B which is natural in the TT-approach; the cutoffs were equipped 
with two different one-parameter famihes of shape functions. 

In 4 dimensions, we found that the Einstein- Hilbert truncation leads to a non-Gaussian 
fixed point for all admissible cutoffs, both of type A and B. The robustness of this predic- 
tion is a nontrivial result since in higher dimensions, for instance, where this truncation 
is less reliable, the fixed point is present for some cutoffs but absent for others. 

Another consistency test successfully passed by the truncation is that all cutoffs agree 
on positive values of and A*. A negative would probably be unacceptable for stability 
reasons, but there is no mechanism in the flow equation which would exclude it on general 
grounds. In fact, < is realized for d < 2. 

For all cutoffs of type A and B the non-Gaussian fixed point is found to be UV 
attractive in both directions of the X-g plane. The linearized flow in its vicinity is always 
characterized by a pair of complex conjugate critical exponents leading to spiral-type 
trajectories which hit the fixed point for k — > oo. This is precisely the stability property 
needed for asymptotic safety. 

By definition, universal quantities are scheme-, or cutoff- independent in an exact cal- 
culation. Truncations lead to a scheme dependence, however. We can use the degree of 
the scheme dependence as a measure for the reliability of the truncation. The critical 
exponents and, as we argued, the product (7*A* are universal quantities. The existence 
of fixed points is a universal feature of the RG flow, but not their precise location in 
parameter space. 

The critical exponents were indeed found to be reasonably constant for a wide range of 
shape parameters. The universality properties of g^^A* are much more impressive though. 
While we found a considerable scheme dependence of g^ and A^, separately, their product 
is scheme independent at a rather amazing level of accuracy. 

We believe that these results hardly can be a mathematical accident, and we consider 
them a very nontrivial confirmation of the hypothesis that the region of parameter space 
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where the non-Gaussian fixed point is situated is well described by the Einstein-Hilbert 
truncation. As a consequence, the fixed point should exist in the exact theory, too. 

Apart from the renormalizability issue the nontrivial fixed point is also intriguing from 



a "phenomenological" point of view. Its relevance for the structure of black holes [T^ and 
the cosmology of the Planck era ||T^ has been pointed out already. Moreover, we saw 
in section |V| that the RG improvement of the graviton propagator suggests a kind of 
dimensional reduction from 4 to 2 dimensions when spacetime is probed at sub-Planckian 
length scales. 



Acknowledgement: We would like to thank C. Wetterich for many helpful discussions. 
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APPENDIX A: THE TT-DECOMPOSITION 



1. Pseudo-projectors for the TT-decomposition 

In subsection [11 B| we introduced the TT-decomposition 



h^.u = + h^, + hi:; (Al) 

valid for arbitrary symmetric rank 2 tensors defined on either closed or open, asymptoti- 
cally flat Riemannian (i-spaces. Here /i^^, /i^^ and /i^^ represent the transverse traceless, 
longitudinal traceless and pure trace part, respectively, which are mutually orthogonal. 
According to eq. ( p.9| ) these parts may be expressed in terms of pure spin-2, spin-1 and 
spin-0 component fields /i^^,, a, (p. In this subsection we show that the component fields 
can be obtained by applying certain operators 11 to the full field hfj_,y. In the following 
these operators will be termed pseudo-projectors. 

As a first step we express the longitudinal-transverse and the pure trace part as 



2 

1 

d' 



hi: = --g,uct^, (P^r^hi:. (A2) 



Here the operator L maps vectors onto longitudinal traceless tensors. Given a tensor hi,, 
the equation {Le)^i, = hi, is solved by 

e^ = i^ + ^D^a (A3) 

where is a transverse vector and a a scalar. This solution is unique up to the addition 
of conformal Killing vectors (CKV's), as discussed in subsection |1I B| . We recover eq. 
by inserting eqs. ( |A2| ) , ( |A3| ) into eq. (^T|). Contrary to e^, the scalar is uniquely 



determined by hfj^^. 

Now taking the covariant divergence of eq. ( [Xl] ) with eq. inserted, and using 

the transversality requirement D^hl, = leads to 
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iVe), = -D^ [h,, - -^g,,cp) (A4) 
with the operator V defined by 

{Ve)^ = -D'iLe)^, . (A5) 
As shown in ref. D is a positive definite, Hermitian operator mapping vectors onto 



vectors. Moreover, the equation {T>s)^ = with an arbitrary given vector always 
possesses solutions which are unique up to CKV's. However, even if these CKV's exist 
they cause no problems in solving eq. ( |A4| ) for e^, see ref. |^D| for details. In order to 



determine this solution we have to invert T). For this purpose we assume that T) (and 
any additional operator that needs to be inverted in the course of this discussion) has a 
complete set of orthogonal eigenf unctions and that the corresponding eigenvalues do not 
have zero as an accumulation point. Then exists and the solution to eq. ( [A^ ) may 
be written as 

e, = {V-^5h), . (A6) 
Here the operator 5 maps tensors onto vectors according to 

= -D^K.^'^^K^ (A7) 

with 

A.."" = \ + 5? - \9,.r' (A8) 
being the operator that projects symmetric tensors onto their traceless part: 

Ku'^^Kp = V - ^9^^u(|) ■ (A9) 



Taking the covariant divergence of the solution ( |A6|) , inserting eq. (|A3| ) , and using -D^^^ = 
yields 

D% = ^D^a = D^{V-'6h)^ . (AlO) 
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This leads to the final result for and a in terms of h^^: 
a = 2{D'^)-^D>'{V-^6h)^ = Eh , 

J, = {V-'6h), - D,{D^)-^D''{V-^~5h), = (Qh)^ . (All) 

Hence the pseudo-projectors 11 which map h^u onto the individual component fields are 
obtained as 

hi, = {UTTh)f,, = (Ah)^, - {LV-^6h)^, , 



^^ = {ULTh)^= i^J-D^-Ricnh 



a = ULLh = y (^2)2 + _A_D^R^^^D, Eh , 

<P = IlTrh = r'h,, . (A12) 

Here we defined 11^7- and Ull in terms of the redefined fields and a which are related 
to and a by eq. (|2]2l|). 

Furthermore, the pseudo-projectors for the transverse decomposition of arbitrary vec- 
tor fields can be inferred from eq. QATTD. For = C^>' + D^{-D^)-2ri with D^.C^'' = 



they are determined by 



c^^ = {nTCf = 0"- D^'i-D^y^D^c . (A13) 

Obviously 11^^ maps tensors onto vectors, Ull and Ilrr map tensors onto scalars, 
and Ul maps vectors onto scalars. Hence these operators cannot be projection operators 
in the usual sense. However, projection operators P mapping arbitrary h^j^^, onto /i^^ or 
h^^, or arbitrary e:^ onto their longitudinal component can be constructed from them 



2|. Since the H's map vectors and symmetric tensors onto their component fields they 



generate a kind of projection in a wider sense of the word. Therefore we call the II's 
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pseudo-projectors. Contrary to 11^^, Ull, ^Tr and 11^, the operators 11^^ and 11^ are 
genuine projection operators mapping symmetric tensors and vectors onto their ST"^ and 
T component, respectively. 



2. Construction of the cutoff and the source terms 

In the present paper we need a formulation for which allows for a description in 
terms of the fundamental as well as the component fields. The translation between the 
two descriptions can be achieved by using the pseudo-projectors for the construction of 
the cutoff and for an appropriate decomposition of the source terms. 

Starting from the definition of the cutoff in terms of the fundamental fields, eq. ( p.25|) , 
we choose the cutoff operators TZf^^ and TZf^ as 

^f= E E nt J7^,)^^,, n,, . (am) 

Here \ (TZk),!,,-)^) - m represents a block of the more general matrix op- 

erator ((^fc)^i^2)^^ ^^^^^^t _ ^} which (7^fc)^i^2 = if both ipi G {C^,r]} and 
■ip2 £ {C'^,v}y or if both ipi G {C^,f]} and 4'2 ^ {C^,v}- The operators (7?.fc),ji^2 
and {'JZk)i(>iip2 are required to satisfy the Hermiticity conditions {T^k)(^2Ci ~ '^'^^)ciC2 
(^fc)^2t/'i = -(^fc)ii^2- Furthermore, we set U^t = Utt, = ^lt, = ^ll, 
= Htv, = Hct = n^, Ilj^ = = n^, and n^, denote the Hermitian 

conjugates of 11^, 11^. Hence the operators 11 appearing in eq. ( |A14| ) can be inferred 
from eqs. ( |AT^) , (IATsI) . 



Inserting eq. (|A14|) into eq. (|2.25|) leads to 



AkS[h, C,C-g] = ^ i^Cih, (7^,)c,c2 ^c. h) 

Cl,(2e{hT,i,a,4>} 



- J2 E {n-,^c,ink)^^,^u^,c) (ai5) 

^l£{CT,fj} i?2e{CT,r;} 
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which can be rewritten in terms of the component fields. In fact, applying eqs. (|A12|) and 
( |M3| ) and using the relation between the matrix operators {T^k)^^^^^ "^^^ i'^k)^-^^^^ stated 
above, we end up with eq. ( p.26|) . 

Let us now consider the source terms in eq. (|2.27|) . Decomposing the sources for the 
ghost fields according to 



K^" = K^,r + D^{-D^)-^^Kfi , K^" = + D^^-D^y-^ (A16) 



with Df,K^r = Df,K^r = yields 



(A17) 



The decomposition of (J, h) is more involved. In analogy with /i^j^ in eq. ( |A1| ) we 
decompose the source J^u into its orthogonal parts: 



J = + + J^' 



(A18) 



According to eqs. (|A2D, ( [A^ ) we write J^^ as 



(A19) 



with 



J^J = (LG)^, = D^e, + D,e, , J^^l; = \{LDv),, = D^D^v - -J^.D^v . (A20) 



1 



Analogous to eq. (|A^), the vector 6^ + D^v /2 with -D^jO^ = is unique up to CKV's. 
According to eq. ( p.l4|) the inner product (J, h) may now be decomposed as follows: 



{J,h) =(/^>■n + (J^/^^n + (J^/^^^ 



(A21) 



It is important to note that, for arbitrary nonvanishing h^^, /i^J, /i^^, /i^^, we have 



(A22) 
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This implies that the operators 11^^, 11^^, Ull, are invertible if their action is 
restricted to fields of the type /i^^, /i^J, /i^^, /i^^, respectively. Therefore the inner 
product of eq. ( |A21|) may be written as 



(J, h) = (j^, u^'tUtt h^) + ( J^ ni^T^LT h"^^) + ( J^ ndULL h 



E (-^cC) • (A23) 



Here we used eq. ( |A12D and introduced the source component 



s 



Jl = ((n^^)^ J^Y = -Ric^^ + 2 (-^2 _ Ric)"^ {D.R'"') 



V 



2 



J<^^(n?,i)tj^^' = .i^'^Vj;. (A24) 



Then combining the results of eqs. (|A17|) and ( |A23| ) we eventually arrive at eq. ( p.28|) . 



APPENDIX B: EVALUATING THE RHS OF THE TRUNCATED FLOW 

EQUATION 

In this section we present several rather lengthy calculations needed for the discussion 
of the Einstein-Hilbert truncation in section |V[ In the following, all calculations are 
performed with g^^ = g^^ where g^y is assumed to correspond to a spherical background 
and the bars are omitted from the metric, the curvature and the operators. 
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1. Computation of the inverse operators (t^^^ + T^k 



' and f + TZk) ' 



In section ^ we derived explicit expressions for the kinetic operators Tj^ 



(2) 



.(2) 



+ 



TZk and S. 



= + TZk- They may be represented as matrix differential operators 



acting on the column vectors {W,^,(l)o,a,(j)i)^ and {v'^ , v'^ , g, g)^ , respectively. In this 
representation they take the form 



/ 



f(2) 
J- u 



[9,9] 



{r'i\9,9] 




hTRT 





02X1 





{^?h9,9] 


O2XI 






^^k\9,9] 



>2xl 



and 



SS![9^9] 



{SS![9^9])^, 




(4? [f'f^ 




V 














(•^gh [^,^] 



0lx2 ^ 
0lx2 
0lx2 

Qk ^ 





(-^gh [^>^] 




(Bl) 



:b2) 



ee 



where 



Qk 



{^k\9.9])^.^ {ft\9,9] 



(B3) 



The entries of these matrices are given in eq. ( |4.16| ). On the RHS of the flow equation ( p.3| ) 
these operators appear in terms of their inverses, which are determined in the following. At 
this point it is important to note that, because of the maximally symmetric background, 
all covariant derivatives contained in the operators ( Pl| ) and ( P^ ) appear as covariant 
Laplacians and that the various entries are x-independent otherwise. This implies that 
these entries are commuting differential operators which allows for particularly simple 
manipulations. Therefore it is not difficult to verify that the inverse operators assume the 
form 
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and 

(•S£'[9.9l)"' 
^ 




with 



{^^"\9,9] 






02x1 





{^^\9,9] 


02x1 






^?\9,9] 

02x1 



0lx2 
0lx2 







1x2 



{SS[9,9])^,^, 




















{sS![9, 



88. 



(•^gh [^>^] 





88 



-1 



frb.j]),,(ff|j.j|),,,.-(rf|j.j| 



(pier 



-[^^^9,9] 



{^^\9,9] 



Inserting these expressions into the RHS of the flow equation (|373|) leads to 



+ -Tr' 
2 



CC. 



CC 



{m9,9])^^{m9,9%^^^-{m9,9]: 

2{r'-^\g,g]+n,[g]) d,{n,[g])^^ 



-Tr' 



E 



1 -1 
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where we used the relations 



(2) 
gh 



vy 



QQ 



o(2) 



q(2) 



1 



5 



5S, 



gh 



uy 



ee 



15 1 5Sffh 



:b81 



The trace of the 0o-term appearing in eq. ( P7|) may be easily evaluated since only the 
scalar eigenmodes and T^'" contribute. We obtain 



It.' 

2 



^ -1 



dt{nk[g]) 



1 A WO) 



2Z 



J2 E d''x^g{x)T'"'ix) Pk + AsMa)R + BsMa)\ 



xdt 



T^'^ix) 



E 



Z^fcA;2i?(o)(A,(d,0)/A;2 



(B9) 



2Zm^ fr^o Az(rf, 0) + Fi?{o)(A;(rf, 0)/P) + Asi(c^, a)R + «)Afc 

Here Ai{d,0) is the eigenvalue with respect to —D^ corresponding to T'"^. Inserting also 
the remaining operators given in eq. ( |4.16| ) into eq. ( [B7D finally leads to eq. ( |4.17| ). 



2. Evaluation of the traces 

In this part of the appendix we evaluate the traces appearing in eq. ( [4.20| ) by applying 
the asymptotic heat kernel expansion. In its original form it has often been used to 
compute traces of operators acting on unconstrained fields, see e.g. refs. For our 

purposes we need the corresponding expansions for operators acting on constrained fields, 
i.e. fields satisfying appropriate transversality conditions. In appendix |^ these expansions 
are derived in detail for Laplacians acting on symmetric transverse traceless tensors, 
on transverse vectors and on scalars, with the following results: 

/ ■ \ "^^^ r 1 
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(25T2) 



^(is-e)D2 
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(rf + l)(rf + 2)(rf-5 + 35rf,2) 



Tr 



(IT) 



-(is-e)D2 



12(d- 1) 



(BIO) 



^47r(s + ie) J 
{d + 2){d-3) + 66d,2 



(fx ^{d-l 



6d 



[is - e)R + 0{R^ 



Tr 



(0) 



,-(is-e)D2 



d/2 



d'^x {1 - - (is - e)R + 0(R^] 
' 6 



(Bll) 



(B12) 



\A7r{s + ie) J 

The next step is to consider an arbitrary function W{z) with a Fourier transform 
W{s). For such functions W, we may express the trace of the operator W{—D'^) that 
results from replacing the argument of W with —D^ in terms of W{s): 



Tr 



W{-D^) = Urn / dsW{s)TT 



-{is-e)D^ 



(B13) 



We obtain the asymptotic expansion of Tr[iy (— D^)] by inserting the heat kernel expansion 
for Tr[e~''^''~^^'^^] into eq. (P13|) . For Laplacians acting on the constrained fields considered 
here they read as follows: 



Tr 



(25T2) 



W{-D')] =(4n)-''/'S^^id-2){d+l)Q,/2[W] j d^x ^ 

+- .^/, —Qd/2^i[W] d'^Xy/gR 



12{d-l] 



Tr 



(IT) 



W{-D')] ={A7Tr''/'Ud-l)Q,/2[W] jd^x^ 



(B14) 



{d+2){d-d,) + Q5d,2 
6d 



Qd/2^i[W] / d'^x^R + Oi- 



, (B15) 



Tr 



(0) 



+lQd/2-i[W] J d^x^R+Oii 



<d~2\ 



(B16) 
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Here the set of functionals Q„[Vr] is defined as 

oo 

Qn[W] = \im [ ds{-is + eyW{s) . (B17) 

e\0 J 
— oo 

By virtue of the MeUin transformation we may now reexpress Qn in terms of W so that 

1 °^ 

Qn[W] = [ dzz^-' Wiz) , n>0; 

1 (n) J 





oo 



Q^[W] = / dzz"'^''-^ ^ ' , n<0, m>-n, mG^ arbitrary. (B18) 

^ ' 

In particular we obtain Qo[W^] = W^(0). 

Let us now consider the case where isolated eigenvalues have to be excluded from 
Tr[W^(— D^)]. According to appendix ^ such traces can be expressed as the differ- 
ence between the complete trace Tr[W^(— D^)] and a term of the form Yl, Di{d,s) 

ie{ii,...,i„} 

xW (A/((i, s)). Here /i, . . . , refer to the modes to be omitted and Ai{d, s) and Di{d, s) 
denote the corresponding eigenvalues of —D^ and their degrees of degeneracy, respectively. 
Since Ai{d, s) oc R we may view W {Ai{d, s)) as a function of R. As outlined in section 
rV| such a function contributes to the evolution of Z^k and only for d = 2, with the 
contribution given by Vr(0). Using the explicit expressions for Di{d,s) (see table 1 in 
appendix 0) and applying eq. (^?7|) we therefore obtain for the traces relevant to the flow 
equation 

Tr'(,^) m-D')] = Tr(i^) [W{-D')] - ^ W{0) J d'x ^R + 0{r<'-') , 
TT'l,^[Wi-D')] = Tr(o)[Py(-D^)] - If W^(O) / d\ ^R + 0{r<'-') , 
Tr'(o) [W{-D')] = Tr(o) [W{-D')] -^-§W{0) J d\ ^R + 0{r<''-') , (B19) 



where the primes have to be interpreted as in section |IV. 

The next step is to insert the expansions of the traces into Sk{R), eq. ( |4.2CI| ), and to 
compare the coefficients of the operators / d'^x^/g and / d'^x^/gR with those on the LHS, 
eq. (|4.5| ). This leads to the following differential equations: 
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+dQd/2 {Ao - 2aAfe)"V - 2dQd/2 [^Wo] | 



(B20) 



t^Nk 



-{2K')'\A7:)-''^L{d) Qdi2-i [(a 



2\i 



-1 



+C2{d) Qd/2-l 

+C4{d,a)Qd/2 



An - 2a\, 



An — 2aXi 



c^{d) Qd/2 



-2 



- 2c2(rf) Q 



d/2-l 



Aq — 2Afc 

-1 



+CM Qd/2 Uo'M) 



35, 



d.2 



_2ZNk{k^ - 2Afc) 2Z^fc(P - 2aXk) 
Here the coefficients q are defined as in eq. 



(B21) 



In eqs. ( P20|) , ( P21|) , the various Qn may now be expressed in terms of the cutoff- 
dependent threshold functions $^ and <l>^ introduced in eq. ( [4.23| ). Using the relations 

Qn [(A + cyM] = A;2("-P+i)$^(c/A;2) - ir/^(A;) p("-f+i)$P(c/A;2) 

Qn [(A + cY^^r,] = k^^--^+'^^l{c/k') (B22) 

we arrive at the differential equations ( |4.21| ) and ( [4. 221 ). 

At this stage the following point should be mentioned. In order to achieve that the 
integrals in eq. ( pi7|) actually converge we have to demand that i?*^''-' {y) rapidly decreases 
as ?/ ±oo. However, since from now on its form for ?/ < does not play a role any 
more we identify {y) with its part for nonnegative arguments and assume that R^^^ (y) 
is a smooth function defined only for y > and endowed with the properties stated in 
subsection IlllB. 



APPENDIX C: FLOW EQUATIONS IN FOUR DIMENSIONS 



In this section we compare our results to those derived in ref. for the d = A 
dimensional case. Inserting d = 4 into eq. ( |4.17| ) we obtain for the RHS Sk of the 
evolution equation 
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^k{R) = 2Tr(25T2) 



^Tr" 



dtPk 



dtPk 



Pk + ^R-2aXk 



dtPk 



A(4,0) 

- 



Pk + - 2«Afc 

dtPkiH^M 



2 



-Tr 



P.(AK4,0))-3^A,J 
Pk + D^ 



Pk — 2Afc 
-Tr 



^Nk 



(25T2) 



t(lT) 



2 -^^(1^) 



Pk 



4 J 



Tr'. 



(0) 



dtPk 

p R 



Pk + D' 



^ 1=0 



4a \ 



- — Tr" 



P,(AK4,0)) 3^_, 
a — 1 



Pk-2\ 



Pk + 



x\ (1 -3«) 



(3 - a)Pk + - — -R 



R - 2aAA. 



+ 4«(« + l)Afe (Pk + D 



-3{l-afJPk\ Pk 



R 



R 



R 



'Pk\lPk-'-^-V^^-D'~ 3 



(CI) 



with Pfc(A,(4,0)) = Ai(4,0) + Pp(o)(Az(4, 0)/P). Our result Q agrees with eq. (3.22) 
of ref. up to a few (typographical) errors occuring in eq. (3.22). To be more precise, 
the prime at the Tr'^j^ji-j-term in lines 1 and 4 of eq. (PT|), the factor Di{A, 0) appearing in 
the first term of lines 3 and 5, and the factor 1/2 contained in the term oc P in line 7 are 



mistakenly left out in the corresponding equation of ref. |21 



Expanding the flow equations for ZjsfkXk and Z]\fk, eqs. ( |4.21| ) and ( [4.22| ), with respect 
to Xk, using the relation dtXk = Zj^\dt{Z]sikXk) + XkrjN{k) and setting, as in |l2ll| . 



= 2k Z 



Nk 



Zu = z 



Nk 



Vk ^ dt In Zk ^ -mik) , = 2$^(0) , eg ^ ^1(0) (C2) 



leads to 



dt^k = (47r)-2p \ l^gl + g + « ) g2 + r^, 



1-1 ^25 .2 



+ OCXk 



(C3) 
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dtXk 



■17 



1 1 

2^2 +^^2^2 



+ 



24 



•47 



13 1 
"24^^ 



24 



(C4) 



In this form the flow equations for the couphngs are suitable for a comparison with the 
corresponding results in eqs. (4. 6)- (4. 9) of ref. Apart from the contributions from 

matter fields, which are also considered there, the results differ by a factor 60/24 in front 
of the gl'terms. Presumably, this deviation can be explained by a wrong sign introduced 
in for a certain term. This term is a contribution from the ghosts produced by heat 
kernel expanding the last two terms in line 3 of eq. ( |U1| ), and it carries the prefactor 
30/24. We may conclude that, apart from these corrections, our results agree with those 
in pl|, which is as it should be since the same cutoff is used. 



APPENDIX D: TENSOR SPHERICAL HARMONICS 

In this section we introduce the spherical harmonics Tj^, Tjf^ and T'*" for symmetric 
transverse traceless (ST^) tensors /i^^, transverse (T) vectors and scalars on a c?- 
dimensional spherical background S'^. These harmonics form complete sets of orthogonal 
eigenfunctions with respect to the covariant Laplacians acting on ST"^ tensors, T vectors 
and scalars, i.e. they satisfy 

-D'Ti::{x) = Md,2)Ti::{x), 

-D'Tl^{x)=Hd,l)Tl^{x), 

-D^ T'"^(x) = Ai{d, 0) r'"(x) (Dl) 
and, after proper normalization, 

5'^ 5- = jd'^x^g {\2ST-)Y'"'" TU^^T^: = jd'^xV^g (I(it))''' T'J^T^'' 

= J (i'^Xv^r'"^T^" . (D2) 
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Here (1(25t2))''''^'' = (d - 2)/2d{g''Pg'''' + gf^'g^P) and (I(it))^'' = {d - l)/dg'"' are the 
unit matrices in the spaces of ST"^ tensors and transverse vectors, respectively. In eq. 
( [Pl| ) the Ai{d, s)'s denote the eigenvalues with respect to —D^ where s refers to the spin 
of the field under consideration and / takes the values s, s + 1, s + 2, ■ ■ -. Furthermore, the 
second upper index at Tj^, Tjj^ and T'™, m, takes the degeneracy of the eigenvalues into 
account. It assumes values from one to Di{d,s) with Di{d,s) the degree of degeneracy. 
In ref. explicit expressions for Ai{d, s) and Di{d, s) are derived which can be found in 
table 1. The eigenvalues are expressed in terms of the curvature scalar R = d{d — l)/r^ 
of the sphere with radius r. In ref. |^B| it is also shown that the spherical harmonics T^™, 
Tjf^ and T'™ span the spaces of ST"^ tensors, T vectors and scalars so that we may expand 
arbitrary functions /i^^, and (p according to 



oo 



A(rf,2) 



1=2 m=l 
oo Di{d,l) 



Im I 



1=1 m=l 



X] 



oo Di{dfl) 

m = T. E 0^mT'-(x). (D3) 

i=0 m=l 

Here the coefficients {/i^}, {iim} and {(pim} are countably infinite sets of constants that 
are uniquely determined by /i^j,, and 0. Eq. ( |D3| ) may now be used to expand also 
any symmetric non-T^ tensor and nontransverse vector in terms of spherical harmonics 
since they may be expressed in terms of ST"^ tensors, T vectors and scalars by using the 
decompositions (|2.9|), ( |2.18| ), see e.g. [P^-|55[]. 

In this context it is important to note that the Di{d,l) = d{d + l)/2 modes T^'™ 
and the Di{d,0) = d + 1 modes T^'™ satisfy the Killing equation (|2.11|) and the scalar 
equation ( p^.l2D , respectively, and that r°'^ = const. As discussed in subsection [11 lj| , 
arbitrary symmetric rank 2 tensors receive no contribution from these modes. In the case 
of arbitrary vectors it is the constant scalar mode that does not contribute. Such modes 
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have no physical meaning and have to be omitted therefore. 



Table 1: Eigenvalues of —D"^ and their degeneracy on the (i-sphere 


Eigenfunction 


Spin s 


Eigenvalue Ai{d,s) 


Degeneracy Di{d, s) 






2 


l{l+d-l)-2 p 


{d+l){d-2){l+d){l-l){2l+d-l){l+d-3y. 


/ = 2,3,... 


d{d-l) 


2{d~iy.{i+iy. 




1 


l{l+d-l)-l p 


l{l+d-l){2l+d-l){l+d-3y. 


/ = 1,2,... 


d(d-l) ^ 


{d-2y{i+iy 







l{l+d-l) p 
d(d~l) 


{2l+d-l){l+d-2y 

i\{d-iy 


/ = 0,1,... 



APPENDIX E: HEAT KERNEL COEFFICIENTS FOR DIFFERENTIALLY 

CONSTRAINED FIELDS 



In this part of the appendix we supply the tools necessary for the evaluation of func- 
tional traces and we derive the heat kernel expansions for Laplacians acting on differen- 
tially constrained fields. 

As a first step we consider a functional trace of the form 



Tr(,[c])[/(-/^^)] 



f{-D 



X 

Ml-. -Ms 

2 \ Mi...Aisi'l...i's 



\ ^ i J'/ Ui...Ustll...tJ.s 



(El) 



x=y 



Here / is an arbitrary smooth function whose argument is replaced with the covariant 
Laplacian defined on the space of spin-s fields with a possible symmetry and/or transver- 
sality constraint C, as indicated by the subscript (s[C]) at the trace. Note that / inherits 
the matrix structure from the corresponding Laplacian. Furthermore, l(s[c]) denotes the 
unit matrix in the space of independent field components. Given a closed^ Riemannian 



WXhe restriction to closed Riemannian manifolds is done for the sake of notational simplicity 
only. In principle, the results extend to noncompact asymptotically flat Riemannian manifolds. 
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manifold {A4,g) we now assume that {f/^^ fisi^)} ^ complete set of orthonormal func- 
tions on {Ai,g) spanning the space of fields under consideration. Then by making use of 
the completeness relation 

d'^ix - y) 1 



eq. (|ET| ) can be written as 



Vl-.-Vs 



(E3) 



Clearly if {V^^ ^^{x)} is taken to be a complete set of orthonormal eigenfunctions with 
respect to the covariant Laplacian such that —D^ U^i...^S-'^) ~ ^^i...^is(^) ^'i- (|E3|) boils 
down to 

Tr(.[c])[/(-I^^)] = E/(A.) • (E4) 

k 

In general the evaluation of such traces is a formidable task and one has to resort to 
approximations. The most familiar such approximation is the "early time" expansion of 



the diagonal heat kernel (x 



x) for |t| —>■ 0. (We write t = s+ie with Im(t) = e > 0.] 



It has been discussed in many references 

d 



X 



X 



^ (4^) ' " ~ + ^ } ■ ^^^^ 

Here Q is an arbitrary smooth matrix potential and the a^'s are tensor polynomials 
proportional to the n/2-th power of the curvature and endowed with the same matrix 
structure as + Q. They depend on the space of fields under consideration. For 
operators + Q acting on unconstrained fields, i.e. fields with independent components, 
the first three coefficients take the form^ 



^§In the context of the Einstein-Hilbert truncation performed in the present paper we only need 
the coefficients ao oc and 02 oc r~^. However, for truncations containing invariants quadratic 
in the curvature we shall also need 04 oc r~^, see |24|. 
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ao(x; Q) = ao = 1 , 
a2{x;Q) = P , 

a^ix; Q) = ^ {R.^a^R^"'"' - R.^R^" + D'r) 1 + + ^T^^uTZ^" + ^^'^ • (E6) 
Here 1 is the unit matrix in the space of field components, and 

P = Q + ].R1. (E7) 
Furthermore, TZ^u is the curvature operator defined as the commutator 

n,,, = [D^, D,] = D^D, - D,D^ . (E8) 
Inserting the asymptotic expansion (^) into eq. (|ETD we obtain for the functional 



trace of the heat kernel 



Tr 



j d'^x ,/g {tr ao - it tr 02(2;; Q) - tr 04 (x; Q) + O (t^) } (E9) 



where tr is the matrix trace with respect to the tensor or spinor indices. 

Let us have a closer look at the coefficients a^. For arbitrary scalars (p, TZp^p vanishes 
since D^D^cf) = D^D^cf). In the case of arbitrary spin-s fields F^^,,,^Xx) with integer spin 
s > 1 we obtain 

TZa/B F^i...fMsi^) = '^Ra/3^,,^ Ff,^...fli_^^lf,i+^...^isi^) ■ (ElO) 
i=l 

Hence, for scalars, 04 receives no contribution from TZajsTZ"'^, while for fields with nonzero 
spin it amounts to a nonvanishing contribution with a {d^ x (i^)-matrix structure. For 
arbitrary vectors and rank-2 tensors we find, respectively. 



fit/ 



(jZapTl"^^ ^^^^ — —Ral3-ytiR°'^^ p gua — Rafi-yu R"''^^ a 9 tip + '^RafinpR'^^ va ■ (EH) 

From now on we restrict our considerations to matrix potentials Q of the form Q = 
qRl(s) with q a real constant, and we assume that the metric corresponds to a maximally 



86 



symmetric background S'^. Setting a„(x; = a„(g) we obtain aQ^q) = oq = l{s) and 

(^2{(l) = (1 + 6q')/6i?l(s), independently of the spin s of the field. Here the dependence on 
s is totally encrypted in the unit matrix This is not the case for the coefficient 04. 
It is given by 



04(9) 



1 
360 



5£ -7d + 6 
d{d- 1) 



+ 60q + 180g^ 



1(0) R 



(E12) 



[«4(g)] 



360 



d5(I3T^ + 60, + 180,' 



(E13) 



[«4(g)] 



fiupa 



1 

360 



h 60g + ISOg^ 



d^{d-l] 



R' 



fiupcr 



3d^{d 



i9^j,u9pu 9p<j9up) 



R' 



(E14) 



for scalars, vectors and rank-2 tensors, respectively, with l(o) = 1, (l(i)j = 9pu and 

(1(2)) = 9^lp9va■ 

Up to this point we considered only unconstrained fields. For fields subject to con- 
straints like transverse vectors and ST'^ tensors /i^j, we cannot directly apply eqs. ([E6| ) 
or ( [E13|) , (|E14|) . However, the heat kernel coefficients for and /i^^, can be computed 
from those of the unconstrained fields using the decompositions ( |2.18| ) , ( |2.9|) for arbitrary 
vectors and arbitrary symmetric tensors h^J^y. 

From appendix we can infer that the sets of orthonormal — Z)^-eigenfunctions 



{T^'"|mG{l,...,AKl)}, / = 1,2,...} 

U { (A;K 0))-^ D;r'^\ m e {i, . . . , aK o)}, / = i, 2, . . .} 



(E15) 



and 



{T^™|mG{l,...,A(^^,2)}, / = 2,3,...} 



U { (2 (ak^, 1) - f 



n t'™ + D T 



Ira 



me {!,..., AKl)}, / = 2,3,... 
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U [Md,0) 



me{l,...,Di{d,0)}, 



/ = 2,3,...|u|^^7^.T'-|mG {1,...,/^, 



'i{d,0)}, 1 = 0,1,... 



(E16) 



span the spaces of all vectors and of all symmetric tensors /i^j^, respectively. Here the 
T'"^'s are the normahzed spherical harmonics of eq. (P2|) . 

Now we insert these eigenfunctions into the trace formula ( |E3| ) with / taken to be 
an exponential. Then we use the commutation relations of appendix |G| in order to pull 
the -D^'s from the TT-decomposition through the exponentials and to combine them 
to Laplacians. This leads to the following decomposition for the traced heat kernels of 
the unconstrained vectors and symmetric tensors in terms of the heat kernels for the 
differentially constrained fields: 



Tr 



(1) 



Tr 



(IT) 



(0) 



(E17) 



Tr 



(2S) 



-it(D2+gR) 



Tr(25T2) 

+Tr(o) 



+ Tr 



(IT) 



-q)R) 



Tr 



(0) 



~it(D2+gi?) 



■-+i)R 



g ^*ld{<i-i) 



hglR 



(E18) 



The last term of eq. ( [E17| ) and the last three terms of eq. ( [E18| ) arise from those 
spherical harmonics T'"^ and T^™ which are not contained in the sets of eigenfunctions 
( [E15| ) and (|E16|) . To be more precise, the last term in eq. (|E17|) comes from the constant 
eigenmode T^'^ of the operator + {dq + l)/dR. Furthermore, the last but second and 
the last but first term in eq. ( [El(^ ) take account of the eigenmodes T°'^ = const and T^'*" 
of the operator + {2/{d — 1) + q)R, respectively. As discussed in subsection [II B| the 
T^'^'s satisfy the scalar equation (|2.12|) and are therefore in a one-to-one correspondence 



with the PCKV's of S . The last term in eq. (|E18D comes from the eigenmodes T^'™" of 



the operator D"^ + {{d + l)/{d{d - 1)) + q)R, which are the KV's of S"^. 
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These subtraction terms compensate for the corresponding unphysical contributions 
contained in the complete traces for the constrained fields on the RHS of eqs. ( |E17| ), 



(E18). This can be seen as follows. Consider the functional trace of eq. (|E3|) . Omitting 
the contributions from the modes Uk^, . . . , t/fc„j we denote the functional trace involving 
only the remaining modes with Tr(gj'(^]-)[/(— Z^^)]. Then eq. ( [E^ ) implies the following 
relation between Tr('j;j'(^]-)[/(— D^)] and the complete tTa.ce Tr(s[q)[/(— Z^^)]: 

Tr'(;['c])[/(-^')] = Tr(,p])[/(-/^^)]- E / p.) • (E19) 

ke{ki,...,k„} 



This rule indeed yields the last term in eq. (|E17| ) and the last three terms in eq. (|E18|) . 

As the next step we insert the asymptotic expansion ([E9|) into both sides of eqs. (|E17|) 
and ( |E18| ) and compare the coefficients of R. This leads to the following Seeley coefficients 
for the constrained fields: 



trCtol(i'jn) 

tr a2(g)|(ir) 



tr flol 



(IT) 



tr flol 



(0) 



tr 04 (g) I 



(IT) 



tra2(g)|(i) - traa 
tr 04(^)1^-,^^ — tr 04 




+ ^6d,2{l + 2q) R^ + ^S,,^R\ (E20) 



tr 



(2ST2) 



tr aol 



(25) 



tr ao 



(IT) 



tra2(g)|(25T2) = 02(^)1 (25) - tra2 



- 2 tr ao 

' d+1 
d{d-l) 



(0) 



+ q 



tr 02 



(IT) 



d-1 



+ q 



(0) 



tr«2(g)|(o) + ^Sd,2R 



tra4(g)|(2ST2) = tr«4(g)|(25) - tra4 



d+1 
d{d-l 

7 



+ q 



tr 04 



(IT) 



d-1 



+ q 



(0) 



tra4(g)Lo) + '5d,2 4 + -g R' + -6d,iR' . 



(E21) 



1(0) ■ \- ■ 2 V ^ '3 
The terms proportional to the 5's originate from the subtraction terms on the RHS of 
eqs. ( |E17|) , ( [E18|) which are due to the unphysical eigenmodes. These terms have an 
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expansion of the form &2m?" while the terms of the heat kernel expansion are of 

m=0 

the form / d'^x^/g tr an oc r'^"". Comparing powers of i? oc only under the condition 

—2m = d — n a given term &2m'"~^'" contributes to / d'^Xy/gti an- Hence for n, m fixed, the 
Seeley coefficients tr a„ for the differentially constrained fields receive a contribution from 
a term of the form 62m '""^"^Z (/ d'^x^/g) at most for one specific value of the dimensionality 
d. In particular, the subtraction terms in eqs. ( |E17|) , ([E18|) do not contribute to troo, 
while tra2 and tra4 on the LHS of eqs. ( |E20| ) , ( IE21D receive contributions from terms of 
the form 6d,2bo/r^ and 5^,2^2/'"'', ^dAbo/r'^, respectively. 

The matrix traces on the RHS of eqs. ([E20|) , (P21| ) can now be evaluated by using 
the heat kernel coefficients for the (differentially) unconstrained fields. For scalars we 



have tra„(g) = an{q) (1 



(0) 



1). For vectors the traces are evaluated according to 



tra„(g) = g^'^[an{q)]fiv so that we obtain from eqs. (^6|), (pill 3D 



tr do If^i'j — d 



tr a2(g)|(i 



l + 6g 



(1) 



6 



dR 



tr 04(9) 1(1 



360 



5d^ - 7rf2 + 6rf - 60 , , 2' 

— + 60dq + ISOdq^ 



(E22) 



d{d-l) 

In order to determine tr a„ for symmetric tensor fields we have to symmetrize the heat 
kernel for unconstrained rank-2 tensors according to 



X 



-islD'^+qR) 



x) 



+ (x 



[D^+qR) 



(2S) 4 

\ fiucrp 



X 



\ pupa 

x) 



(2) 



+ X 



^-is(D^+qR) 



Xj 



+ {X 



(2) 



[D^+qR) 



\ upcrp 
X ) 



(2) 



\ vppa 

x) 



(2) 

(E23) 



before we can apply eq. (|E9| ) with eqs. (|E6|) and ( [E14|) . This leads to 



[«n(g)] 



pupa 



(25) =i([««(^)] 

+ [an{q)] 



pupa 



(2) 



[an{q)] 



uppa 



puap 



(2) 



+ Mq)] 



upap 



(2) 
(2) 



(E24) 
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and in particular {l(2S))fiupa = {gfipdua + gfiadup) /2. For tensors the matrix traces are 
computed according to tra„ = g'^''g'^'^[an{(l)]fiupcT which yields 



tr ao|(-25') 
tr a2(g)|(2s) 



-d{d+l) 



1 + 6g 
12 



d{d+l)R 



tra4(g)|(2s) = ^ 



720 



5d^ - 2d^ - d'^ - lUd - 2A0 
d{d- 1) 



+ 60q{l + 3q)d{d+l] 



Finally we insert the matrix traces eqs. ([E22|) , ( [E25| ) and tri 



1(0) 



I (0) 



(E25) 



into eqs. 



21|) and determine the heat kernel coefficients for transverse vectors and ST"^ 



tensors. The results are summarized in table 2. 

Let us add a final remark concerning the applicability of the asymptotic expansion 
([E5|). Since it is valid only in the limit |t| ^ it is clear that it cannot be integrated over 
Re(t) = s or Im(t) = e term by term, in general. However, this is possible if the heat 
kernel is integrated against a "test" function which suppresses large values of s or e. This 
is indeed the case for our application of the asymptotic expansion presented in subsection 



B2 



Table 2: Heat Kernel Coefficients 


field 


tr flo 


tr a2{q) 


tr 04 (g) 


tensor 


{d-2)id+l) 
2 


(d+l)(d+2)(d-5+3 5d,2) p 


(d+l)(5(i*-22(i3-83d2-392(i-228+1440(5d,2+3240 5d,4) j-,2 


12{d-l) 


720(i(d-l)2 

, ^ (d+l)(d+2)(d-5+3 5d,2) d2 , ^2 (rf-2)(d+l) d2 




I2(d-l) \ q ^ n 


T 
vector 


d-1 


(d+2)(d-3)+6<5d,2 D 


5d'^~12d^-47d'^-186d+180+360 5a,2+™5a,4 r,2 


+q{d-l)R 


360d2(d-l) ^ 
, „ (d+2){d-3)+6<5d,2 d2 , „2 d-1 d2 

+ q R +q ^ R 


scalar 


1 





5ci2_7d+6 p2 1 1^ 02 1 1^2 p2 
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APPENDIX F: VARIOUS COEFFICIENTS 



In this appendix we define the coefficient functions which appear in eqs. ( |4.13| )- (|4.2CI|) 
of subsection p.V C] and in eq. ( |B5| ) of appendix [B 1| . 

, , d{d~3)+A , , a(rf-2)-l , , a{d - A) 



d{d- 1] 



d 



. /, ^ a{d-2)-2 ^ , d-A , , 

As2{d,a) = — — — — — , Assid) = —— , As4{d,a^ 



2a{d-l) - (rf-2) 
_a{d-2)-2 



As5 = 



a{d - 2) - 2{d - 1) ' ' d 

{d -2){d + 2)c? + [d^ - lOrf + 8)a + 2[d - 2) 



d 



d'^a 



Bsi{d,a) = 



2ad 



2a{d-l) -{d-2) 



Bs2{d,a) = 



2ad 



2aid-l)-id-2) d-2 
Csi{d,a) = -— — — — — — - — -, Cs2{d,a) 



a{d - 2) - 2{d - 1) ' 

d- 12(d- 1) -a(c/-2) 



4(d- 1) -2a{d-2)d- 1 



rf2 



a 



Cs3(,d,a) = 

Fsi{d,a) = 
Fs2{d,a) = 



Es{d,a) = -^^{2a{d-l)-{d-2)) 



a{d - 2) - 2{d - 1) ' ' ' 4«rf2 

2 [a{d - 2) - 2{d - 1)] [2a{d - 1) - (rf - 2)] 

4(rf-2)(c/-l)(a-l)2 



d^ 



a 



(Fl) 



APPENDIX G: COMMUTATION RELATIONS FOR A MAXIMALLY 

SYMMETRIC BACKGROUND 

In the following we summarize the commutation relations which were used in order to 
derive eq. ( |4.13| ) of subsection [IV C| , and eqs. (|E17| ) and ( [E18| ) of appendix ^ They are 



valid for the class of maximally symmetric backgrounds. 

D,D,e = ^iu (Gl) 
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(D^t + D,i^) [D'^C + = (^D^ + f) + divergence (G3) 

1 ^o_\ /^„^„_ 1 



= — ; — aD^ \ + I a + cov. divergence (G6) 

a \ a — 1 J 

{d,1 + D,i,) exp [d^) {d'^C + D'^e) 

= -2i, I^D' + exp I^D' + e + cov. divergence (G7) 

D^a^ exp (f)^) D^'a = -a exp + ) ^ + divergence (G8) 



exp(D") ( D^D'^ - i^^'^D" ) a 



d-l^-.f-n R \ f 2R \ ^ , , 

— —aD^ + exp + a + cov. divergence (G9) 

a \ a — 1 \ a — 1 



APPENDIX H: APPROXIMATE SOLUTIONS FOR THE FIXED POINT 

In the following we determine the approximate formula for the position of the non- 
Gaussian fixed point discussed in subsection |V C| . In a first approximation we set = 



A,. = and determine g^, from the condition r^^r* = 2 — d alone. Solving this equation for 
g^, leads to 

^* Bi{K-a,d)-{d-2)B2{K]a,d) ^ ' 

which, for A* = 0, boils down to 
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2-d 



9* 



(47r)f-i<| h{d) $^/2_i(0) + hid) <fi/2-i(0) + hid, a) 



d/2- 



■d/2\ 



+hid,a) $2/2(0) + 35^,2 



.l-2Afe l-2aAfc 

Here h, ■ ■ ■ are d- and a-dependent coefficients defined as 

,,,, rf^ - 4rf3 - 9^2 _ 12 id - 2)id^ - 13d^ - 2Ad + 12) 
hid) = TTTTTl — r\ ' ^2(«) = 



(H2) 



hid, a) 
hid, a) 



12rf(rf- 1) 
d^ - 4d^ + 9d^~8d~2 



2Adid - 1) 



- id~2)a 



2rf(d- 1) 

id-2)id^ -Ad^ + 5d'^ -8d + 2) (rf-2)2 



-a . 



(H3) 



4rf(rf - 1) 2 

Employing the exponential shape function (|3.8| ) with s = 1, and setting d = 4 and a = 1, 
for instance, eq. (|H2|) yields (7* ~ 0.590. Here we used that for this shape function 
$1(0) = TrVe, $^(0) = 1, $1(0) = 1, $2(0) = 1/2. 

As a different approximation scheme, we determine (A*,(7^,) from a set of Taylor- 
expanded /3-functions. Using 

d 



dXk 
d 



$^(-2«Afc) =2«p$r^(-2aA 



dX 



$^(-2aAfc) = 2ap$^+i(-2«A, 



(H4) 



we expand the /3-functions (|4.28|) and ( [4.30|) about gk = Xk = and obtain 



(3^iXk, gk] a, d) = -2Xk + i^ddgk+ kdid - 1 + 2a)(47r2)^-^ "^IM - id - 2) cuJ A^ gk 



+ldid+l)id - 2)iAn)'-iu,^l/,iO)gl + O {g') 



f3JXk, gk] a, d) = id- 2) gt-id- 2)ujd gl + O (g 



(H5) 



Here z/^ and are defined as in eqs. ( ^?B| ) and ( p.ll| ), and O (g^) stands for terms of 
third and higher orders in the couplings gi(/c) = A^ and g2(fc) = gk- Now g^: is obtained 
as the nontrivial solution to f3g = 0, which reads 
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9* = = ^(4vr)^^ {hid) $^/2_i(0) + hid, a) <^1/M}~' ■ (H6) 

Inserting eq. (|H6|) into /3;^ and neglecting also the terms quadratic in the couplings the 
condition [3^ = leads to 

A* = ^ = $1/^(0) {k,id) $i/,_,(0) + hid, a) $^/,(0)}-^ . (H7) 

Using the shape function (|3.8|) with s = 1 we obtain from eqs. (|H6|) and ( p7| ) in c? = 4 



dimensions 



_ /ISvr 55 



which yields iX*,g*) = (0.287,0.751) for a = 1, for instance. 
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